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This report describes a measurement of the top quark mass in pp 
collisions at a center of mass energy of 1.8 TeV. The data sample was col- 
lected with the CDF detector during the 1992-95 collider run at the Fermilab 
Tevatron, and corresponds to an integrated luminosity of 106 pb _1 . Candi- 
date ti events in the "lepton+jets" decay channel provide our most precise 
measurement of the top quark mass. For each event a top mass is determined 
by using energy and momentum constraints on the production of the tt pair 
and its subsequent decay. A likelihood fit to the distribution of reconstructed 
masses in the data sample gives a top mass in the lepton+jets channel of 
176.1 ± 5.1(stat.) ± 5.3(syst.) GeV/c 2 . Combining this result with measure- 
ments from the "all-hadronic" and "dilepton" decay topologies yields a top 
mass of 176.1 ±6.6 GeV/c 2 . 
PACS numbers 14.65.Ha, 13.85.Qk, 13.85.Ni 



Section 1 



INTRODUCTION 

This paper describes a measurement of the top quark mass using 
events produced in proton-antiproton (pp) collisions at the Fermilab Tevatron 
with a center-of-mass energy of 1.8 TeV and reconstructed through the de- 
cay mode ti —>■ W + h + W~b —>■ i + vh + qq'h (and charge conjugate mode). 
Throughout this paper the symbol I will be used to denote either an electron 
or a muon exclusively. We present results from two data samples with inte- 
grated luminosities of 19.7 pb -1 (Run la) and 86.3 pb _1 (Run lb) collected 
with the Collider Detector at Fermilab (CDF) from September 1992 to June 
1993 and from February 1994 to July 1995, respectively. 

The existence of the top quark was established by direct experimental 
observation at the Fermilab Tevatron by the CDF Jl], ^ and D0 collabora- 
tions ||]. These analyses led to ti cross section and top quark mass mea- 
surements. Additional analyses showed that the kinematics of the observed 
events were inconsistent with being solely from background sources and were 
consistent with standard model ti j|. With substantially larger data samples 
and improved understanding of systematic uncertainties, more precise mea- 
surements of the top quark mass || Q and ti production cross section |7|, || in 
pp collisions were recently reported. The larger data samples were used to per- 



1 



2 



form detailed comparisons of kinematic variables between ti candidate events 
and simulated standard model ti and background events ||. The data sam- 
ples were also used in the identification and analysis of ti production into fully 



hadronic final states [10, 11] and final states involving two leptons, 11 [L2L O] 



or £t 18 



The top quark is defined as the I3 = +1/2 member of a weak SU(2) 
isodoublet that also contains the b quark. In pp collisions, top quarks are 
expected to be produced primarily in ti pairs via quark-antiquark annihilation 
(~90%) or gluon fusion («10%) and decay through the electroweak interaction 
to a final state consisting of a W boson and b quark. In the standard model, the 
branching fraction for t — > Wb is expected to be nearly 100%. The decay width 
is calculated to be 1.6 — 1.7 GeV for masses between 150 and 180 GeV/c 2 [14|. 
The top quark mass is sufficiently large that top-flavored hadrons are not 



expected to form [15 



The mass of the top quark, M top , is an important parameter in cal- 
culations of electroweak processes since it is approximately 35 times larger 
than that of the next heaviest fermion. Like other fermion masses, M top is 



not predicted in the standard model Ul6| . On the other hand, the standard 
model relates the masses of the top quark and W boson to that of the Higgs 
boson, so that precise measurements of the former imply bounds on the lat- 
ter. With the assumption of the validity of the standard model, experimental 
studies of the electroweak interaction can alternatively be used to estimate the 
value of M top . For instance, a fit to LEP (including LEP-II) data, leaving the 
top quark mass and the Higgs mass as free parameters, yields an inferred top 
quark mass of 160^g 3 GeV/c 2 and a Higgs mass of 60±lf GeV/c 2 (TJ. 

The decay modes of the W bosons into either lepton-neutrino (£u), 
{tu) or quark-antiquark (qq') final states classify candidate ti events into four 
main categories. All-hadronic final states, which comprise approximately 44% 
of ti decays, correspond to those events in which both W bosons decay hadron- 
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ically. Lepton+jet events are those events in which only one of the two W 
bosons decays hadronically while the other decays into iv and form 30% of 
ti decays. Dilepton events are defined as those in which the W bosons decay 
into either ev or fiu final states and occur only about 5% of the time. Lastly, 
there is an additional 21% of events for which the final state includes one or 
more r leptons. The r events are particularly difficult to identify because r's 
decay into leptons or hadrons and are often indistinguishable from the other 
final states, thus contaminating the other samples. Each ti decay mode is 
characterized by a final state consisting of two b hadrons and either zero, two, 
or four additional jets, depending on the decay mode of the W's in the event. 
Additional jets beyond those from the ti decay may also arise from initial and 
final state radiation of the incoming and outgoing partons. 

The direct experimental determination of M top through analysis of 
ti pairs produced in pp collisions can be obtained by comparing observed 
kinematic features of top events to those predicted for different top quark 
masses |J. While any kinematic variable which exhibits sensitivity to the 
mass of the top quark may be used to measure M top , the lowest statistical 
uncertainty is achieved by explicitly reconstructing the top mass from the ti 
daughter decay products. In this paper, we discuss the complete reconstruc- 
tion of top events in the lepton+jets topology and report the measurement 
of M top obtained using the distribution of the reconstructed top quark masses 
from the data sample. 

This paper is structured as follows. Section ^| presents a descrip- 
tion of the CDF detector, emphasizing the subsystems most important to this 
analysis. Section [I] discusses the reconstruction of jets and leptons in the CDF 
detector and defines the sample of events which are used in the measurement 
of the top quark mass. Section f| describes the simulations used and discusses 
the details of the background calculation. Section |5| describes the corrections 
which are applied to the raw calorimeter measurements. Section |6] presents 
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the algorithm used to estimate the top quark mass on an event-by-event basis 
and describes the results of the algorithm when applied to simulated samples 
of both tt and background events. The description of the likelihood procedure 
and the subsequent extraction of M top are the subjects of Sections [7| and |[ 
Section | describes the systematic uncertainties associated with the top quark 
mass measurement. Combining the measurements from the lepton+jets, dilep- 



ton, and all-hadronic analyses is the focus of Section [10]. Conclusions are given 
in Section [Ll]. 



Section 2 
THE CDF DETECTOR 



The CDF detector is an azimuthally symmetric general purpose de- 
tector. It consists of independent subsystems designed for distinct tasks. The 
three most relevant subsystems to tt detection are the tracking chambers, the 
calorimetry, and the muon chambers. In this section, we briefly describe these 
subsystems. The various subsystems are shown in the side view of one quad- 
rant of the detector in Fig. [2.1| . A more detailed description of each of these 
components can be found in Refs. |], [L9| 



2.1 Detector subsystems 

The tracking system consists of three subsystems that are all im- 
mersed in a 1.4 T solenoidal magnetic field. The outermost system, the central 
tracking chamber (CTC) 0], is a wire drift chamber consisting of 84 concen- 
tric cylindrical layers of sense wires. The CTC has a length of 3.2 m and an 
outer radius of 1.32 m which results in full acceptance for charged particles 



in the region \r]\ < 1 |21|. The momentum transverse to the beamline (Pt) is 



measured by the CTC with a precision given by <5(Pt)/Pt = O.OOIIPt (Pt 
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Figure 2.1: Side view of one quadrant of the CDF detector for Run 1. The 
detector is symmetric about the interaction point. 
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in GeV/c), when the track is constrained to go through the beam position 
determined for each run. 

Inside the CTC is a set of time projection chambers (TPC) with 
tracking coverage in the region |?7| < 3.25. This detector, referred to as the 
VTX, is used to measure the position of the pp interaction vertex along the 
z-axis with a resolution of 1 mm. In events with more than one reconstructed 
vertex, the primary vertex is taken to be the one with the largest number 
of VTX hits on its associated tracks. Primary collisions are spread with an 
approximately Gaussian density along the z-axis with a ~ 30 cm. The primary 
vertex is required to be within ±60 cm of z — 0.0. The efficiency of this 
requirement is evaluated using the same techniques described in Ref. |2B[ and 
is estimated to be 95.6%. 

The innermost tracking system, the silicon vertex detector, SVX, 
consists of four layers of single-sided silicon detectors (the Run la detector 
was replaced for Run lb due to radiation damage) P4"| , mounted inside two 
cylindrical barrels having a combined length of 51.0 cm. The four layers are 
located at radii of approximately 3.0, 4.2, 6.8 and 7.9 cm from the beam- 
line. The axial strips of the three innermost layers have 60 /im pitch, and 
the outermost layer has 55 /im pitch. The silicon detector measures hits in 
the transverse plane with a precision of 13 /im and the impact parameter of 
tracks relative to the primary vertex has a precision of (13 + 40/Pt) /im (Pt 
in GeV/c). Secondary vertices (from weak decays, for example) are identified 
and reconstructed by augmenting reconstructed CTC tracks with hits found 
in the SVX. The precision of the SVX enables efficient identification of sec- 
ondary vertices from the decays of b hadrons (cr ~ 400/iin). The momentum 
resolution of a track reconstructed using both the SVX and CTC detectors is 



given by <5P T /P T = ^/(0.0009P T ) 2 + (0.0066) 2 , where P T is in GeV/c and the 
second term is due to multiple scattering. 

Muons are identified by the association of reconstructed track seg- 
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ments in the proportional wire chambers of either of the three muon systems, 
the central muon system (CMU) [|25|| , the central muon upgrade (CMP), or 
the central muon extension (CMX), with charged particle tracks observed in 
the CTC The CMU and CMP, separated by 0.6 m of steel, each cover the 
pseudorapidity region \r]\ < 0.6. In that region CMU covers ~ 84% of the solid 
angle, CMP « 63%, and both combined ~ 53%. At larger pseudorapidities 
the CMX provides ~ 71% coverage of the solid angle for 0.6 < < 1. 

The CDF calorimeters are segmented into projective towers. The 
towers are further divided into compartments designed to separately measure 
electromagnetic and hadronic energy. Three separate regions of calorimetry 
provide coverage in r\ from —4.2 to 4.2. All of the electromagnetic calorime- 
ters use lead as the absorber, while the hadronic calorimeters use iron. In the 
central region, coverage is provided by electromagnetic and hadronic calorime- 
ters, CEM |26| and CHA/WHA |7|, respectively. The CEM is composed of 
alternating layers of lead and scintillator, whereas the CHA and WHA are 
composed of alternating layers of iron and scintillator. Coverage at larger pseu- 
dorapidities is provided by calorimeters PEM and PHA, and in the far forward 
regions by the FEM and FHA. These calorimeters use gas proportional cham- 
bers, instead of scintillators, as the active sampling medium. The calorimeters 
provide identification of, and energy measurement for jets, electrons, photons, 
unclustered energy ||28|| , and missing transverse energy (Et ) @. The cover- 
age in pseudorapidity and the energy resolution for the calorimeters are given 
in Table O. 



2.2 Luminosity and triggers 

The events used in this analysis are extracted from two data samples 
with integrated luminosities of 19.7 pb~ 1 (Run la) and 86.3 pb _1 (Run lb) 
collected during the period from September 1992 to June 1993, and from 
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Table 2.1: Coverage in pseudorapidity and energy resolution for the vari- 
ous calorimeters. The symbol © signifies that the constant term is added in 
quadrature with the sampling (first) term. Energy resolutions for the electro- 
magnetic calorimeters are for incident electrons and photons. For the hadronic 
calorimeters, they are for incident pions. Et should be expressed in GeV. 

February 1994 to July 1995 respectively. Instantaneous luminosities varied 
between 1 x 10 30 to 2 x 10 31 cm 2 sec 1 during the data taking period, with 
averages that increased from w 3.3 x 10 30 cm~ 2 sec~ 1 during Run la to ~ 
1 x 10 31 cm~ 2 s~ 1 for Run lb. The corresponding average number of interactions 
per crossing increased from 0.6 to 1.8. Since the measured jet energies increase 
in the presence of additional interactions, the corrections to the jet energies 
differ between Run la and Run lb (see Section |5.1.1|) . 

A multi-level trigger is used to select events containing high-Px lep- 
tons P], [/J. To increase the tt acceptance in the muon channel, a trigger based 
on the missing transverse energy (Et ) was added for Run lb 0. For the 
high-Px inclusive lepton sample, only triggers from the central region are used 



in this analysis. The CEM trigger efficiency for fiducial [|30| electrons from tt 
events with Et > 20 GeV and \q\ < 1 is essentially 100%. The muon trigger 
is measured to be 85.4% efficient for fiducial muons from tt events that have 
P T > 20 GeV/c. 



Section 3 
DATA SAMPLES 



The data sample selection for this analysis is based on standard model 
decay of top quark pairs through the tt — > iuqq'bbX channel. The final state 
should therefore include a high-E T (Pt) electron (muon), significant missing 
transverse energy and four jets. The momenta of these objects are measured 
from data recorded with each detector subsystem, sometimes in combination. 
The four-momenta of electrons are expressed in terms of (Ex, <f), r], m) where 
Et is the transverse energy (Et = Esin9), <p is the azimuthal angle, 77 is the 
pseudorapidity and m is the mass. For muons and jets Pt is used rather than 
Et- In all cases, the direction of these objects is measured with much greater 
precision than their energies. In this section, we first describe the identification 
and reconstruction of leptons and jets, and then we define the data samples. 



3.1 High-Px leptons 

We are most interested in identifying charged leptons which are pro- 
duced from the decay of a W boson. These leptons are distinguished from 
those produced in semileptonic decay of b or c quarks because leptons from 
W-boson decay are not part of a jet and have typically much higher Pt- A 
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sample of high-Px leptons is used to select leptons which are consistent with 
having come from H^-boson decay. 

A sample of events which contain high-Ex electrons are selected from 
the Run 1 data sample by requiring the electron to have Et > 20 GeV/c and 
be in the central region of the detector (|^| < 1). Backgrounds from photon 
conversions and charged hadrons are rejected by cutting on several variables. 
Here we describe those cuts which provide the largest discrimination against 
background. A detailed discussion of other selection criteria can be found in 
Ref. 0. Electrons are required to have a CTC track pointing to the electron 
shower in the CEM. The energy in the hadronic calorimeter divided by the 
energy detected in the electromagnetic calorimeter (HAD /EM) is required 
to be less than 5%. We also require that the energy of the shower divided 
by the momentum of the associated track is less than 1.8. Electrons from 
photon conversions are removed with an efficiency of 88% [Q] by requiring 
each electron to have a matching track in the VTX, and the invariant mass 
of this track with any other CTC track to be greater than 0.5 GeV/c 2 . The 
energy of high-Ex electrons is measured using the calorimeter energy in the 
tower to which the CTC track points plus the adjacent towers High-Ex 
electrons are measured with a resolution of <r(Ex)/Ex = 13.5%/a/Ex © 2%, 
where Ex is in GeV. 

The high-Px muon sample is selected by requiring that each event 
contain at least one muon candidate which has Px > 20 GeV/c and is in the 
central region of the detector (\rj\ < 1). Muon candidates are identified by a 
match between a track segment in CMU, CMP, or CMX and the CTC. The 
primary backgrounds are from secondary particles in charged hadron show- 
ers which "punch through" the calorimeter and produce tracks in the muon 
chambers, and cosmic rays. To reject the charged hadron background, the 
muon is required to have an energy deposition in the calorimeters which is 
characteristic of a minimum ionizing particle. Backgrounds from cosmic rays 
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are rejected by requiring that the track extrapolates back (in r — <p) to within 
3 mm of the beamline and that in the r — z plane it is within 5 cm (at r=0) of 
the primary vertex. A number of other selection requirements are made which 
are described in Ref. [l|]. The momentum of high-Px muons is measured by 
constraining the CTC track to the average beam position. Its transverse mo- 
mentum is measured with a resolution of ct(Pt)/Pt — 0.11%Pt, where Pt is 
in GeV/c. 

From these high-Px lepton samples, we further select those events in 
which the high-Px lepton is isolated |32| from jet activity. For the lepton+jets 
analysis, we require that there is only one W — > Iv candidate in the event. 
The lepton I is referred to as the primary lepton in the event. 



3.2 Jet reconstruction 

Jets are constructed from calorimeter tower information using a cone 
algorithm with cone radius AR = Ar/ 2 + A<p 2 = 0.4. The jet transverse 
energy is defined as the sum of the energy deposited in calorimeter towers 
within the cone, multiplied by sin#, where 9 is the polar angle of the Ex- 
weighted centroid of the clustered towers. After correcting for the various 
energy losses (see Section ||), jets which do not contain heavy flavor, and have 
Pt > 80 GeV, have a transverse momentum resolution of <5Pt/Pt ~ 12%. A 



discussion of the jet reconstruction algorithm can be found in Ref. [33, 49 . 



3.2.1 Identification of 6-quark jets 

The identification of jets that arise from b quarks (6-quark jets or 
simply b jets) plays an important role in the analysis described in this report. 
The identification relies on finding evidence for a I?-hadron decay, using two 
separate tagging algorithms. 

The silicon vertex (SVX) tag algorithm [|1], [7| searches within a jet for 
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displaced vertices due to 5-hadron decays. It is applied to jets that have raw 
Et > 15 GeV and uses tracks which are within AR < 0.4 of the jet axis and 
have hits in the silicon vertex detector. The algorithm allows for two passes. 
In the first pass, a secondary vertex is required to have at least three tracks 
with Pt > 0.5 GeV/c, at least one of which has Pt > 2.0 GeV/c. In the 
second pass, tighter track quality cuts are applied, and a secondary vertex is 
required to have at least two tracks with P T > 1.0 GeV/c, including at least 
one with Pt > 2.0 GeV/c. This algorithm has an efficiency of about 48% for 
tagging at least one b jet in a tt event. 

The soft lepton tag (SLT) algorithm [[l], [54j searches for additional 
leptons which are consistent with having come from a semileptonic I?-hadron 
decay. The lepton is required to have Pt > 2 GeV/ c and to be within AR < 0.4 
of a jet with raw jet E T > 8 GeV. The efficiency for tagging at least one b jet 
in a tt event with this algorithm is about 15%. 

The SVX algorithm obtains both higher purity and higher efficiency 
then the SLT algorithm. However, the SLT algorithm is also employed for 
tagging b jets because it uses nearly uncorrelated information and adds signif- 
icantly to the acceptance. 

3.3 Top mass candidate sample 

Full reconstruction of candidate tt events is possible if the event has at 
least four jets and a W candidate decaying into either ev or \w. The majority 
of such events are not from tt production but rather from the production of 
a W boson in association with jets. The fraction of these background events 
containing at least one b jet is of the order of 1% Q, while standard model 
tt decays are expected to always have two b jets. Data samples with larger 
fractions of tt events can therefore be formed by requiring evidence of b hadrons 
in one or more jets. 
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To facilitate the measurement of the top quark mass, we apply se- 
lection criteria which are expected to increase the fraction of ti events in the 
sample. We refer to these events as the Top Mass Candidate Sample, and they 
satisfy the following cuts: 

1. High-Ex lepton trigger satisfied; The event should have an electron 
(muon) with E T > 20 GeV (P T > 20 GeV/c) and \r}\ < 1. 

2. Et , as calculated using the raw tower energies, is greater than 20 GeV. 
For events with a primary muon this % T includes a correction for the 
muon momentum. 

3. The candidate primary electron or muon track must be isolated and 



of good quality (see Section |3.1| ). Only one isolated lepton should be 
present. 

4. Candidate dilepton (tt — > V r vt~vbbX) events, defined according to the 
selection criteria of Ref. [12 1, are rejected. 



5. Events with Z-boson candidates are removed. A Z-boson candidate is 
defined by two oppositely charged, same flavor high-Px leptons (Pt > 
20 GeV/c) that have an invariant mass between 75 and 105 GeV/c 2 . 
Also, we remove the event if it includes a high-Px photon |35[ and the 
invariant mass falls in the Z mass window. 



6. The primary vertex of the event must be within 60 cm of z — 0.0. 

7. At least three jets with E T > 15 GeV and \r)\ < 2.0. 

8. For events with exactly three jets satisfying criterion 7 above, we require 
at least one additional jet with E T > 8 GeV and \rj\ < 2.4. 

9. After the mass reconstruction is performed, events are required to pass 
a goodness-of-fit cut, x 2 < 10.0, where the variable \ 2 is defined in 
Section 0. 
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A sample of 324 events pass criteria 1-7, and are the same as those 
used in the CDF measurement of the tt production cross section [|7J. Criteria 1- 
9 are identical to those used in all our previous measurements of the top quark 
mass |l], @- After imposing criteria 1-8, our sample consists of 163 events. The 
last requirement removes 12 events, from which we obtain an inclusive sample 
of 151 V^+multi-jet events. Thirty-four of the events have SVX or SLT tagged 
jets. As discussed below, the Top Mass Candidate Sample is estimated to 
consist of approximately 74% background. Requiring the presence of 6-tagged 



jets improves considerably the signal-to-background ratio (see Section p. 3.1 



3.3.1 Mass subsamples 

To describe the mass subsamples which are used in this analysis, it 
is helpful to decompose the Top Mass Candidate Sample into two classes of 
events which are expected to have different signal-to-background ratios (S/B). 
Class I events have exactly three jets with Et > 15 GeV and \r)\ < 2 and one or 
more additional jets with Et > 8 GeV and \tj\ < 2.4. Class II events have four 
or more jets with Et > 15 GeV and \rj\ < 2. Because of the larger amount of 
energy contained in the four leading (i.e., four highest Et) jets, class II events 
have a larger S/B than class I. 

Previous measurements of the top quark mass at CDF used a com- 
bined sample of 6-tagged events [[T], 0] that contained events from both class 
I and class II. Monte Carlo simulations show that the statistical uncertainty 
on the measured top quark mass is reduced by 10% by combining the results 
of separate fits on three non-overlapping subsamples of events. The first sub- 
sample consists of events that have one and only one SVX tag. The second 
subsample consists of events in which there are two SVX tags. The third one 
includes events that have one or two SLT tags, but no SVX tags. Further 
Monte Carlo studies show that an additional 7% improvement is obtained by 
including the No Tag events from class II. The 75 No Tag events excluded from 
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the Top Mass Candidate Sample are expected to have a background fraction 
of 93%. Inclusion of these events does not improve the statistical uncertainty 
on the top quark mass measurement. To summarize, the four mass subsamples 
are [06| : 



• SVX Double : Events with two SVX tags; 

• SVX Single : Events with one and only one SVX tag; 

• SLT : Events with one or two SLT tags, but no SVX tags; 

• No Tags : > 4 jets with E T > 15 GeV and \r]\ < 2. 

The numbers of data events in each of these subsamples are shown in Table |3.1| . 
In categorizing the events into the subsamples, tags are only counted if they 
are on one of the four highest Et jets. This choice is made because the four 
leading jets are assumed to be the primary partons from the ti decay (see 
Section |6|). Also shown in the table are the expected S/B ratios, using the 
background estimates presented in Section |4.3.3|. The measurement of the top 



quark mass in the lepton+jets channel is based on these four subsamples. 



Data Sample 


Number of Events 


Expected S/B 


SVX Double 


5 


24 


SVX Single 


15 


5.3 


SLT 


14 


0.8 


No Tags 


42 


0.4 



Table 3.1: Subsamples used in the lepton+jets mass analysis and the expected 
signal to background ratio (S/B) for each. See Section [4.3.3| for background 
estimates for these subsamples. 



Section 4 



SIMULATION AND 
BACKGROUND 



This section describes the Monte Carlo methods used to simulate the 
signal and background events, and the estimation of the background in the 
four mass subsamples. For this purpose we use Monte Carlo programs that 
generate the signal and background processes contributing to the data sample, 
and a detector simulation which models the response of the detector to the 
final state particles. Unless otherwise noted, the Monte Carlo programs use 
the MRSDO' |37| set of structure functions. Detailed properties of 6-hadron 



decay, based on observations from the CLEO experiment ||38|| , are included in 
all the Monte Carlo generators. The response of the detector to the final state 
particles is parametrized using distributions observed in data. See Section |5] 
for details on the calorimeter simulation. 



4.1 Signal modeling 

The simulation of ti events relies mainly on the HERWIG [[39| (Ver- 
sion 5.6) Monte Carlo program. Additional checks are provided by both 
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PYTHIA gg (Version 5.7) and IS A JET glj (Version 6.36). HERWIG is based 
on the leading order QCD matrix elements for the hard process, followed by 
coherent parton shower evolution, cluster hadronization, and an underlying 
event model based on data. PYTHIA is similar to HERWIG in that it is based 
on leading order QCD matrix elements; however, partons are fragmented using 
the Lund string model. ISAJET is a parton shower Monte Carlo program based 
on the leading-order QCD matrix elements for the hard-scattering subprocess, 
incoherent gluon emission, and independent fragmentation of the outgoing 
partons. 



4.2 Background modeling 

The Monte Carlo program used to study the kinematics of the back- 
ground is VECBOS This is a parton-level program based on tree- level 
matrix element calculations for VT+jets production. The simulated events 
produced by VECBOS contain a W boson and up to four additional final state 
partons. These partons are subsequently evolved and hadronized using a sep- 
arate program [f43[ derived from the parton shower model contained in the 



HERWIG Monte Carlo generator. The CDF simulation program is then used 
to simulate the detector response and produce the final sample of background 
events for further analysis. 

The VECBOS events generated for this analysis use the W + 3 parton 
matrix elements, with the required additional jet being produced during parton 
showering. The Q 2 scale of the hard scatter is set to the square of the average 
P T ((P T ) 2 ) of the outgoing partons unless otherwise noted. 

The VECBOS Monte Carlo generator has been shown to reproduce 
distributions of a wide range of kinematic variables in a large sample of IV+jets 



events [44] in this experiment. In addition, distributions of kinematic variables 



have been studied in tt-depleted and tt-enriched subsamples of W+ > 3 jet 
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events in this experiment ||. The Monte Carlo simulations reproduce the 
distributions in both subsamples when we use the expected fractions of HER- 
WIG (for tt) and VECBOS (for background) events. Further checks which 
demonstrate that VECBOS is appropriate for background modeling are given 



in Section 3.5 



4.3 Background estimation 

In the measurement of the top quark mass, we constrain the fraction 
of background events in each of the mass subsamples to an expected value. 
The computation of the expected value for each mass subsample is achieved 
by first computing the expected number of background events from relevant 
background processes for both class I and class II events (see Section |3.3.1|) . 
Some of the background processes are computed as absolute predictions while 
others are given as a fraction of the number of background P^-candidates in the 
data sample. The expected tt and background fractions (which sum to unity) 
in the Top Mass Candidate Sample are then estimated by using a maximum 
likelihood fit which compares the observed rates of events with SVX and SLT 
tags with predicted rates. The predicted rates, which use estimates of the 
tagging probabilities for tt and background events, depend on these fractions. 
The tt fraction is a free parameter in the fit, and is allowed to vary to optimize 
the agreement between the observed and predicted numbers of tagged events. 
The fitted tt fraction in the Top Mass Candidate sample is then combined with 
SVX and SLT tagging probabilities to evaluate the expected tt and background 
contribution in each of the mass subsamples. The same principle has been used 
to measure the ti cross section using W+ > 3 jets events J7| (£5"] . 

The tagging probabilities we use, and the contributions of various 



background channels, are similar to those in Ref. [45], but are not identical 



because of differences in the event selections and the exact tagging rules. The 
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event selections used in this paper require a fourth jet and impose a x 2 cut on 
the kinematic mass fit (described in Section []). The tagging rule used here, 
requires that the SVX and SLT tags are counted only if they are on one of the 
four leading jets in the event. The resulting differences in tagging probabilities 
and backgrounds are determined using the HERWIG and VECBOS Monte Carlo 
simulations. 

4.3.1 Inputs into the background calculation 

The inputs into the calculation are the background processes, their 
expected rates, and the corresponding SVX and SLT tagging probabilities. 
The rates and tagging probabilities are estimated for both the class I and 
class II events of the Top Mass Candidate Sample. Of the 151 events in the 
Top Mass Candidate Sample, 87 are in class I and 64 in class II. 

The background processes are classified into two categories: contribu- 
tions which are computed as an absolute number of events, and contributions 
which are calculated as a fraction of the number of background candidate 
iy+jets events (N w ) in the data sample. In the latter case, the contribu- 
tion includes Z+jets events that pass the lepton+jets selection criteria. The 
background processes considered are listed in Tables |4.1| and [4.2| for the two 
classes (the processes are the same for both classes). The expected numbers 
of background events from the different processes are also given in the tables. 

For the first six processes we have absolute predictions. For the 
W+jets and Z+jets processes we have predictions for each process relative to 
their sum. The last two columns in Tables |44] and (L2| give the SVX and SLT 
tagging probabilities per event for each background process. The probabilities 
in rows 1-13 are for cases where there is a real displaced vertex or a real soft 
lepton. Each of the background processes can also contribute fake SVX and 
SLT tags (mistags), and these probabilities are given in row 14. In either case, 
the SVX and SLT tagging probabilities include the requirement that the tag 
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is on one of the four leading jets and take into account the \ 2 cu ^ on the 
kinematic mass fit. 



Item 


Background 


Number of 


tsvx 


e SLT 


# 


process 


events 


(%) 


(%) 


Absolute backgrounds 


1 


non-W/Z 


5.7±0.8 


4.3±2.2 


2.5±1.8 


2 


WW 


0.7±0.1 


5.8±1.7 


1.3±0.7 


3 


wz 


O.liO.O 


5.8±1.7 


1.3±0.7 


4 


zz 


O.OiO.O 


5.8±1.7 


1.3±0.7 


5 


Z — > TT 


0.9±0.1 


3.5±2.5 


4.6±4.6 


6 


Single Top 


0.4±0.1 


30.6±7.0 


9.0±2.4 


W/Z+jets backgrounds 


7 


Wbb 


(0.028 ± 0.004)^ 
(0.056 ± 0.013)iV^ 
(0.053 ± 0.016)^ 
(0.005 ± 0.002)iV^ 
(0.005 ± 0.002)iV^ 
(0.001 ± 0.001)iV^ 
s 0.85N^ 


22.7±3.1 


7.0±1.9 


8 


Wcc 


5.7±1.0 


5.5±1.2 


9 

10 


Wc 
Zbb 


3.7±0.5 
22.7±2.0 


6.3±1.8 
7.0±1.9 


11 


Zee 


5.7±1.0 


5.5±1.2 


12 


Zc 


3.7±0.5 


6.3±1.8 


13 


W/Z + u, d, 


0.0 


0.0 


Mistag probabilities 


14 


1-13 




0.4±0.1 


3.2±0.4 



Table 4.1: Backgrounds which contribute to class I events in the Top Mass 
Candidate Sample. Shown are the contributing processes, their estimated 
contribution, and the SVX and SLT tagging probabilities per event for each 
process. Backgrounds whose absolute rate is calculable (a total of N* bs events) 
are given by 1-6. Backgrounds that are given as fractions of the number of 
W/Z '+jets events in the data sample are given by 7-13. Nfy is the total 
number of W/Z +jets background events in class I. All background processes 
contribute to SVX and SLT mistags, with the probabilities listed in row 14. 
There are 87 events in class I. 

The expected backgrounds and tagging probabilities are calculated 
as follows. The non-W/Z background is calculated directly from the data [0. 
The WW, WZ, and ZZ background rates are evaluated by multiplying the 
acceptances for these processes as determined from the PYTHIA Monte Carlo 
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Item 


Background 


[ Number of 


tsvx 


tSLT 


# 


process 


events 


(%) 


(%) 


Absolute background calculations 


1 


non-W/Z 


5.5±1.7 


4.3±2.2 


2.5±1.8 


2 


WW 


0.7±0.2 


5.8±1.7 


1.3±0.7 


3 


WZ 


O.liO.O 


5.8±1.7 


1.3±0.7 


4 


zz 


O.liO.O 


5.8±1.7 


1.3±0.7 


5 


Z — > TT 


0.7±0.3 


3.5±2.5 


4.6±4.6 


6 


Single Top 


0.3±0.1 


30.6±7.0 


9.0±2.4 


W/Z+ jets backgrounds 


7 


Wbb 


(0.054 ±0.012)iV^ 


27.4±2.7 


7.5±2.6 


8 


Wcc 


(0.087 ±0.025)JV# 


6.0±1.0 


5.6±1.2 


9 


Wc 


(0.073 ±0.022)^ 


3.8±0.5 


6.3±1.8 


10 


Zbb 


(0.003 ± 0.003)iV^ 


27.4±2.7 


7.5±2.6 


11 


Zee 


(0.003 ±0.003)^ 


6.0±1.0 


5.6±1.2 


12 


Zc 


(0.001 ± 0.001)iV# 


3.8±0.5 


6.3±1.8 


13 


W/Z + u, d, 


s 0.78JV# 


0.0 


0.0 


Mistag probabilities 


14 


1-13 




0.4±0.1 


4.2±0.5 



Table 4.2: Backgrounds which contribute to class II events in the Top Mass 
Candidate Sample. Shown are the contributing processes, their estimated 
contribution, and the SVX and SLT tagging probabilities per event for each 
process. Backgrounds whose absolute rate is calculable (a total of N^ s events) 
are given by 1-6. Backgrounds that are given as fractions of the number of 
W/Z+jets events in the data sample are given by 7-13. is the total 

number of W/Z+jets background events in class II. All background processes 
contribute to SVX and SLT mistags with the probabilities shown in row 14. 
There are 64 events in class II. 
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simulation by their production cross sections [KJ. The Z — > tt background 
is estimated using the PYTHIA Monte Carlo simulation. The normalization is 
obtained by scaling the number of reconstructed Z — > ££+ > 1-jet events in the 
simulation to the number observed in the Run 1 data sample. For single top 
quark production, we use the PYTHIA and HERWIG Monte Carlo programs 
to evaluate the acceptances for the W* — > tb and VT-gluon fusion processes 
respectively. The production cross sections are normalized to the published 
theoretical values fETl. 



The expected fractions of Wbb and Wcc events in the data sample 
are evaluated using the HERWIG and VECBOS Monte Carlo programs. For 
each jet multiplicity bin, the expected background is given by the product of 
the corresponding background fraction, tagging probability and the number 
of VT-candidate events. The Wc background is estimated from HERWIG in 
an analogous way to what is done for the Wbb and Wcc backgrounds. The 
Zbb, Zee and Zc backgrounds are calculated using a combination of HERWIG, 
PYTHIA and VECBOS. The simulations show that in both the Z+l jet and 
Z+2 jet multiplicity bins Zbb events are approximately twice as likely to pass 
our kinematic cuts as Wbb. The corresponding ratio for Zee to Wcc is approx- 
imately 1, and Zc/Wc is about 0.3. We assume that these scalings also hold 
in the higher jet multiplicity bins. The Zbb, Zee and Zc background rates are 
thus obtained by scaling the Wbb, Wcc, and Wc rates by 2.0 ± 0.5, 1.0 ± 0.3, 
and 0.3±0.15 respectively. The overall Z/W normalization is determined from 
the data sample, and is 0.092±0.020 for events in class I and 0.030±0.030 for 
events in class II. 

The SVX and SLT tagging probabilities in lines 1-13 in Tables [4.1| 



and [4.2| give the probability per event, that one or more jets will be tagged 
due to the decay of a long-lived particle (i.e., a b hadron, a c hadron, or a 
r). For backgrounds which are computed using Monte Carlo programs, the 
tagging probabilities are evaluated by simulation of the detector's response 
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to the final state particles of each of the background processes. For SVX 
tags, the probabilities are calculated using only jets which have an uncorrected 
E T > 15 GeV and \rj\ < 2. For SLT tags the probabilities include all jets which 
have an uncorrected Et > 8 GeV and \r]\ < 2.4. The tagging probabilities for 
W + u, d, s are set to zero since these events have a negligible contribution 
from long-lived particles. 



The SVX and SLT mistag probabilities (line 14 in Tables |4.1| and 
are estimated by applying "mistag-matrices" to the jets in each event 
of the Top Mass Candidate Sample. The mistag matrices [|lj] for SVX and 
SLT tags are measured from inclusive jet data and describe the probability 
for a jet that does not contain heavy flavor to be tagged by the SVX and 
SLT algorithms respectively. Monte Carlo simulations show a lower mistag 
rate in background events than in tt events, with a ratio of 0.70±0.05 for 
both SVX and SLT tags. This ratio is included in the mistag probabilities 



shown in Tables |4.1|44.3| . The effect of using equal mistag probabilities for tt 
and background has been investigated, and the resulting background numbers 
change by a negligible amount. 

Tagging probabilities for tt events were determined using the HER- 
WIG Monte Carlo program. Additional checks of these probabilities were 
provided by both the PYTHIA and ISAJET simulations. The probabilities for 
tagging at least one 6-quark jet in a tt event are shown in Table [4.3| . Also 
shown are the probabilities for tagging a jet which does not contain heavy 
flavor (mistags). As before, the SVX and SLT tagging probabilities include 
the requirement that the tag is on one of the four leading jets and require the 
X 2 cut on the kinematic mass fit. 

4.3.2 tt and background fractions in each event class 

We first estimate the fractions of background and tt events in each 
of the two event classes defined in the preceding section. For each event class, 
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Tagging probabilities per tt event 




e S vx(%) 


esM%) 




Class I Class II 


Class I Class II 


Real tags 


44.8±4.5 49.9±5.0 


14.9±1.5 14.8±1.5 


Mistags 


0.6±0.1 0.7±0.1 


4.8±0.5 6.4±0.7 



Table 4.3: SVX and SLT tagging probabilities in tt events for class I and 
class II events. Shown are the probabilities for tagging one or more jets which 
contain b or c quarks (real tags) and the probabilities for tagging one or more 
jets which do not contain b or c quarks (mistags). 

we compare the expected rates of tags with the observed rates in each of 
four subsamples. The subsamples are events with (i) only SVX tags, (ii) only 
SLT tags, (iii) both SVX and SLT tags, and (iv) No Tags. The division into 
these subsamples was chosen to optimize, according to Monte Carlo studies, 
the background fraction estimate, and is not identical to the mass subsample 
division. Note that for subsample (iii) the tags can be on the same jet or on 
different jets. 

The expected numbers of events in each of these subsamples (indexed 
by j) can be calculated as a function of the numbers of tt events (iVjf) and 
non-top VT+jets events in the event class, using an expression of the form: 

N; xp = a, x Nu + £ c) x N* bStj + £ d)b) x N w (4.1) 

k i 

Here the first term gives the expected contribution from tt events, and the last 
two terms give the expected number of events from background processes. The 
indices k and i refer to the background processes 1-6 and 7-13, respectively, in 
Tables |4.1| and [4.2| . The parameter aj is the (SVX or SLT) tagging probability 



for ti events in the j th subsample, while Cj and bj are the tagging probabilities, 
including those for mistags, for background processes k and i. The quantities 
represented by dj are the coefficients of N^r and in Tables [iT] and [12 
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The parameter N% bs ■ is the expected number of background events from the 
k th process. Equation |47T] applies separately to both class I and class II events. 

The tagging probabilities in the expression above are derived from 
the values in Tables |4.1| - |4.3| , apart from some correlation terms. Correlation 
terms between real and mistag probabilities and between SVX and SLT tag 
probabilities are included in the calculations, but these terms are relatively 
small and their effect on the final result is negligible. 

To determine the background and tt contributions to class I and class 
II events, we constrain the total number of tt and background events (i.e., 
summed over the subsamples) to be equal to the observed number of events in 
each class. Then we have just one parameter for each class, the fraction, f t t, of 
tt events (or, equivalently, the fraction of background events). A given value 
of ftt determines values of N t i and Nw to be used in Eq. [4TT| . A maximum 
likelihood method is used to determine a best estimate of The likelihood 
has the form: 



£=n *$(/*) ( 4 - 2 ) 

i 

where the i th event falls into subsample j and the expected fraction of events 
in subsample j is F l -{f t t)- 

The results of the maximum likelihood fit are f t i = 0.13j?og(siat.) ± 
.01(syst.) for class I, and f t i = 0A5t[\l(stat.) ± .05(syst.) for class II. The 
statistical uncertainties correspond to changes in In L from the maxima by 0.5 
units. The systematic uncertainties result from adding in quadrature the many 
contributions due to changing all the relevant input rates and probabilities one 
at a time by their stated uncertainties. These f t i values imply that ti events 
comprise 11. of the 87 class I events and 28.5^7^ of the 64 class II events. 



The numbers of tt and background events are summarized in Table 4.4 



To check that the model we are using is reasonable, we compare the 
expected numbers of events in each subsample with the observed numbers. 
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The comparison is presented in Table fO|, and shows reasonable agreement 
between expected and observed numbers. 



Process Class I Class II 



tt 

Absolute Backgrounds 
iy/Z+jets 



11.5+& 28.5±f;l 
7.9 ±0.9 7.4 ±1.8 
67.6l|| 28. lt^ 



Table 4.4: Estimated composition of the Top Mass Candidate Sample for class 
I and class II events using the background likelihood fit described in the text. 
Shown are the expected contributions from tt events, absolute backgrounds 
(as listed in lines 1-6 in Tables 4T and 4.2 ) and W/Z+]ets events. The sum of 
each column is constrained to the number of observed events in the Top Mass 
Candidate Sample. 



Observed Expected 



Subsample 


Class I 


Class II 


Class I 


Class II 


Only SVX tags 


3 


10 


5.6 


12.4 


Only SLT tags 


6 


8 


4.2 


4.8 


Both SVX and SLT tags 


3 


4 


1.1 


3.0 


No tags 


75 


42 


76.0 


43.8 


Total 


87 


64 


87 


64 



Table 4.5: The number of observed and expected events in the four subsamples. 
The expectation values are based on the background likelihood fit described 
in the text. The events are separated into class I and class II events. 



4.3.3 tt and background events in the f+jets mass sub- 
samples 

Having found the numbers of tt and background events for the sam- 
ples in class I and class II, we can go to the next step, i.e., compute the 
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expected numbers of top and backgrounds events in the mass subsamples. To 
arrive at estimated ti fractions in the mass subsamples, we need probabilities 
for two SVX tags in an event. We must also combine the ti fractions for class 
I and class II events in each tagged subsample. The untagged mass subsample 
only contains class II events. 

For most of the background channels the probabilities for two real 
SVX tags (i.e., tags due to 6-hadron, c-hadron, or r decays) are very small or 
zero. The non-negligible probabilities are given in Table |4.6| . Our calculations 
for the SVX Double subsample do allow appropriately for real and fake tags 
in all channels. The probabilities for events to enter into one of the four mass 
subsamples use the probabilities in Tables P~T| , fO| , and [O], and are computed 
as follows: 

P(SVX Single) = P(SVX) - P(SVX Double) (4.3) 
P(SLT (no SVX)) = P(SLT) - P(SVX ® SLT) (4.4) 
P(No Tag) = 1 - P(SVX) - P(SLT (no SVX)) (4.5) 

The symbol ® in the second line is used to signify the probability of obtaining 
both an SVX and SLT tag in the same event. 





Double SVX Tag 


Probability per Event (%) 


Process 


Class I 


Class II 


W{Z)bb 


1.9±0.5 


3.7±1.0 


ti 


12.0±2.4 


16.4±3.2 



Table 4.6: The probability per event to have two SVX-tagged jets for W(Z)bb 
background processes and for ti events. Double SVX-tag probabilities for all 
other background processes are negligible and are set to zero. The probabilities 
are evaluated for class I and class II events. 

The computation of the expected ti fraction in each of the mass 
subsamples proceeds as follows. First, for each mass subsample, we calculate 
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the expected ti fraction in each event class. Then, the tt fractions for class I 
and II events are combined into a single tt fraction. For each class, the expected 
tt fraction, g m , in mass subsample m is given by the following expression: 

Afm 

m _ tt,exp / 4 g x 

NT\ ' 

1 * tot, ex p 

An expression of this form applies to both class I and class II events in each 
mass subsample. The numerator, N$ is the expected number of tt events in 
mass subsample m, and the denominator is the expected total number {tt + 
background) of events. The expected total number of events in subsample m 
is calculated using an expression of the form shown in Eq. |4.1| (replace j with 
m, and use the tagging probabilities appropriate for the mass subsamples). 
The tt fractions for each event class in mass subsample m are then combined 
into a single tt fraction, fff, using the following expression: 

Jtt N? + N% ' 1 ' 1 

Here, N™ and Nf} are the observed numbers of events, and gf and gj} are the 
predicted fractions of tt events in the two event classes in subsample m. The 
expected number of ti events is given simply by the numerator of Eq. [O. For 



the No Tag mass subsample we have / t f = g™j, because only class II events 
contribute. Table |T7| shows the observed number of events, the expected total 
number of events (ti + background) and the expected contribution from ti 
alone. Note that the total number of ti events in each class is the same as that 
of Table [O] as expected. 

The /S 1 have both statistical and systematic uncertainties. The sta- 
tistical uncertainties are asymmetric, and are convoluted with the systematic 
uncertainties separately for classes I and II, and the results are in turn con- 
voluted. The systematic uncertainties are assumed to be Gaussian. The end 
result is a likelihood function for each f£r, which is used in the mass likeli- 
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Class I Class II 





Total 


Total 


tt 


Total 


Total 


tt 


Subsample 


Obs. 


Exp. 


Exp. 


Obs. 


Exp. 


Exp. 


SVX Double 


3 


1.5 


1.4 


2 


4.8 


4.7 


SVX Single 


3 


5.3 


3.7 


12 


10.7 


9.6 


SLT (no SVX) 


6 


4.2 


1.2 


8 


4.8 


2.8 


No Tags 


75 


76.0 


5.2 


42 


43.8 


11.4 


Total 


87 


87 


11.5 


64 


64 


28.5 



Table 4.7: The number of observed events, N b s , in the mass subsamples, the 
total expected number of events, and the expected number of ti events. Events 
in class I with No Tags are not used in the top quark mass analysis. 

hood fit described in Section [7|. These negative-log-likelihood distributions as 



a function of the expected number of background events are shown in Fig. [11 
Finally, the estimated composition of each mass subsample can be 
calculated from the fjg values and the various tagging probabilities and event 



rates. The result is shown in Table 4.8. The contributions from mistags are 



included in the sums for each process. From the Table we see that 80% of the 
background is from iy+jets and Z+jets, and another 15% is from non-W/Z 
events, i.e., from multijets (including hh events). The remaining 5% is from 
diboson events, Z — ► tt, and single-top production. The background fraction 
per subsample varies from 4% for SVX Double tagged events to 73% for No 
Tag events. 
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Item 




svx 


SVX 


SLT 






# 


Process 


Single 


Double 


(no SVX) 


No Tags 


Total 


1 


non-W/Z 


0.5 


0.0 


1.0 


4.6 


6.1 


2 


WW 


0.1 


0.0 


0.1 


0.6 


0.8 


3 


wz 


0.0 


0.0 


0.0 


0.1 


0.1 


4 


zz 


0.0 


0.0 


0.0 


0.1 


0.1 


5 


Z — > TT 


0.1 


0.0 


0.2 


0.5 


0.8 


6 


Single Top 


0.2 


0.0 


0.1 


0.2 


0.4 


7 


Wc+Zc 


0.2 


0.0 


0.8 


1.7 


2.7 


8 


Wbb + Zbb 


0.8 


0.2 


0.4 


1.1 


2.5 


9 


Wee + Zee 


0.4 


0.0 


0.8 


2.0 


3.2 


10 


W/Z + u,d,s 


0.2 


0.0 


4.1 


19.6 


23.9 


Background sum 


2.4i°;? 


0.2 ±0.1 


7.6 ± 1.3 


30.4H? 


40.7 


11 


tt 


12.6 


4.8 


6.4 


11.6 


35.3 


Observed events 


15 


5 


14 


42 


76 



Table 4.8: Expected composition (in events) for the four mass subsamples 
from various processes. The Vr+jets and Z+jets processes have been summed 
together. The No Tag subsample only includes contributions from class II, as 
only these events are used in the top quark mass analysis. 



32 




SLT (no SVX) NoTag 



Figure 4.1: The negative log-likelihood function for obtaining a given number 
of background events in each mass subsample: (a) SVX Double tags, (b) SVX 
Single tags, (c) SLT (no SVX) tags, and (d) No Tag events. 



Section 5 



CORRECTIONS TO RAW 
CALORIMETER ENERGIES 



Calorimeter information is used to estimate the jet momenta and the 
net transverse momentum of the particles recoiling against the tt system. This 
section details how those estimates are made. The signal from each calorime- 



ter tower is converted into a raw |48| energy estimate. Tower energies are then 
used to evaluate the total energy in the event and other quantities used in the 
top mass analysis. The raw measurements are corrected for non-instrumented 
regions, non-linear response of the calorimeter, multiple interactions at high 
luminosity, and other effects, before a constrained fit is applied to the ti can- 
didate events. Also in this section checks of the jet energy scale are discussed, 
this being the source of the largest systematic uncertainty in the measurement 
of the top quark mass. 



5.1 Jet corrections and 

The raw momentum of a jet is 
momenta from all the towers belonging 



their uncertainties 

calculated by adding vectorially the 
to the jet cluster (see Section |3~2| ). 
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Tower momenta are calculated from tower energies with the assumption that 
they are energies of particles with zero mass [fE| that originate from the re- 
constructed primary vertex and are located at the center of the tower. To 
measure the top quark mass from candidate ti events, corrections are applied 
to the raw jet momenta in order to obtain estimates of the momenta of the 
daughter partons in the ti decay. The corrections occur in two stages. 



• A set of "flavor-independent" corrections is applied to all jets with 
raw E T > 8 GeV. 

• A second set of corrections, specific to ti events, is applied to the leading 
four jets which are assumed to be the daughter jets from the ti decay. 
These corrections are applied after the flavor-independent corrections, 
and map the measured jet momenta to the momenta of the partons in 
the ti decay. 

A description of the corrections to the raw jet momenta is the focus of this 
section. 



5.1.1 Flavor-independent jet corrections 

To account for detector and reconstruction effects, raw jet transverse 
momenta are corrected using a set of "flavor- independent" jet corrections |[33| . 
The following expression includes all the corrections applied: 

P T (R) = (P£ aw (R) x f rel - UEM(R)) x f abs (R) - UE(R) + OC(R). (5.1) 

The parameter R=^ (Arj) 2 + (A(f>) 2 is the cone radius chosen for the jet mea- 
surement; R=0.4 for this analysis. The corrections are described below: 

• f re i, the relative energy scale, corrects for non-uniformities in calorimeter 
response as a function of rj. 
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• UEM(R) takes into account energy due to multiple interactions in the 
event. 

• fabs(R), the absolute energy scale, maps the raw jet energy observed 
in a cone of radius R into the average true jet energy. This average is 
determined in the central calorimeter assuming a flat Pt spectrum. 

• UE(R) takes into account the energy due to the underlying event, i.e., 
the energy from the primary pp interaction due to fragmentation of par- 
tons not associated with the hard scattering, 

• OC(R), corrects for the energy expected to be outside the cone radius 
of 0.4. 



The f a bs{R) and the OC(R) corrections are functions of the transverse mo- 
mentum of the jet. The relative correction is primarily dependent on the 
pseudorapidity of the jet, with only a weak dependence on the jet momentum. 

The reconstruction of jets starts with the raw clustered energy, P™ w (0.4). 
An uncertainty of ±1% is assigned to the stability of the calorimeter over the 
course of the data taking period. This systematic uncertainty was evaluated 
by comparing the response of the calorimeter to single charged tracks between 
data from Run 1 and data from the 1988-1989 run, which was used for the en- 
ergy calibration discussed later. No systematic difference was observed. Also 



the raw inclusive jet cross section |50| obtained with the 1988-1989 data run 
was compared with that of the Run la data (after correcting for multiple in- 
teractions) and it was found that the ratio was consistent with unity at the 
5% level. Because of the rapidly falling Et spectrum, this corresponds to an 
upper limit on a difference in the energy scale of 1%. 

The relative correction is derived from di-jet balancing data and cor- 
rects for the relative response of the different calorimeter sections to that 
of the calorimeter in the central region (0.2 < \rj\ < 0.7) [59]. The plug 
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Uncertainty on 
77 1 range relative correction 



0.0 - 


0.1 


2.0 % 


0.1 - 


1.0 


0.2 % 


1.0- 


1.4 


4.0 % 


1.4- 


2.2 


0.2 % 


2.2 - 


2.6 


4.0 % 


2.6- 


3.4 


0.2 % 



Table 5.1: The percentage uncertainty on the relative jet energy correction 
for various detector rj ranges. The cracks in 77 between different detectors are 
located near 77=0, 1.2, and 2.4, and have larger uncertainties than the regions 
away from the cracks. 



(1.1 < < 2.4) and forward (2.4 < < 4.2) regions are thus calibrated. 
The precision to which this calibration is known is limited mostly by the num- 
ber of di-jet events available. The effects of different resolutions of the central 
and plug calorimeters on the energy measurements were studied using Monte 
Carlo simulation and are properly included. The uncertainty is larger near 
the cracks between the different detectors due to smaller statistics and worse 



energy resolution. Table |TT] gives the uncertainty (in %) on the relative cor- 
rections for various detector 77 ranges. 

The corrections for multiple interactions (UEM) in the same event 
and the underlying event (UE) in the primary interaction are derived from 
minimum bias data. The average number of interactions in Run la (N„=0.6) 
is different from that of Run lb (N„=1.8), hence a different procedure is used 
for the two samples. For the Run la sample, 0.72 GeV/c is subtracted from the 
jet Pt after the absolute correction and accounts for both effects on average. 
For Run lb, the effects of the underlying event and additional interactions 
are separated. To account for multiple interactions, prior to the absolute 
correction, 0.297 GeV/c is subtracted from the jet Pt for each additional 
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reconstructed vertex in the event. This correction is obtained by studying the 
amount of energy in the event as a function of the number of vertices over the 
course of the run. For the underlying event (UE), we subtract 0.65 GeV/c 
from each jet after the absolute correction. 



The uncertainty on the UEM correction is estimated to be 100 
MeV/c for each vertex in the event. The uncertainty in the UE correction 
is evaluated by looking at variations in the energy density at ±90° with re- 
spect to the two jets in di-jet events when varying the maximal Ex threshold 
on the third jet from 5 to 15 GeV. Based on these studies, we assign a ±30% 



relative uncertainty to the underlying event correction [[49], [33]. For jets with 
Pt > 20 GeV/c the uncertainty is typically less than 0.5% of the jet's Pt, as 
shown in Fig. |5.1| . 



The absolute correction is derived from data and Monte Carlo plus de- 
tector simulation. The simulation includes many features of the CDF calorime- 
ters, the main ones being: non-linearity, cracks and less sensitive regions, sin- 
gle tower thresholds. The response of the calorimeter to incident pions and 
electrons is studied using testbeam data, minimum bias runs, special runs 
which triggered on events containing single isolated tracks, as well as stan- 
dard data runs. The detector simulation has been tuned to agree with these 
data. The step from individual particle response to jets is achieved by tuning 
the Monte Carlo (ISAJET) fragmentation parameters to reproduce a number 
of distributions observed in di-jet data: number of charged particles, spectra 
and invariant mass of charged particles, and the ratio of charged to neutral 
energy f3"3]| . The derived correction then accounts for non-linearity of the 
calorimeter, energy losses near the boundaries of different calorimeter wedges, 
response variation as a function of the position along the wedge and all the 
other effects included in the simulation. The absolute correction, / a 6 s (0.4), as 
a function of corrected jet Pt, P^ r , is shown in Fig. [5.2| (a). 
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Figure 5.1: Uncertainty in jet Ex scale as measured with a jet clustering cone 
of size 0.4. The vertical axis shows the extent to which the measured jet Ex 
response varies due to different systematic effects. 
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Figure 5.2: "Flavor-independent" jet corrections, for a jet clustering cone of 
R=0.4. (a) Absolute correction, f abs , and out-of-cone correction factor, 1+ 
OC/Pt, versus corrected jet P T , Pf"" • 0) Total correction, P™ r /P r T aw , as a 
function of Pf r - (c) total correction, P™ r /P™ w , as a function of P™ w . (d) 
fraction of measured momentum, P?j? u '/P™ r versus PJj? r . 
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The systematic uncertainty in the absolute correction is attributed 
to (a) calorimeter response, and (b) fragmentation related effects 



sec 



Fig. fj.ip . The parameters that describe the calorimeter's response to incident 
electrons, photons and pions have uncertainties due to finite statistics and 
assumptions which are made. For example, at low momentum (\p\ < 5 GeV/c), 
the largest source of uncertainty in the charged pion response comes from 
the estimation of the amount of energy in the shower from 7r°'s. Additional 
uncertainty comes from the uncertainty in the relative response across the face 
of a calorimeter cell and the energy deposition in cracks between calorimeter 
cells. The uncertainty in the calorimeter's response to photons is assigned to 
be the same as for electrons. Uncertainty in the fragmentation parameters 
comes from the modeling of the tracking efficiency in jets, and the level of 
agreement between the simulation and data. 

The contributions to the jet Et uncertainty from these sources are 
evaluated by shifting the input values of these parameters by +1 and — 1 
standard deviation (+lcr and —la), and calculating the resulting shift in the 
reconstructed jet energies. For (a) we separately vary the pion, electron, and 
photon responses by +la and —la, and add the resulting shifts in the jet 
energies in quadrature. For (b), we vary the charged tracking efficiency by its 
uncertainty and reevaluate a new set of fragmentation parameters. These new 
fragmentation parameters are in turn varied one at a time, and the resulting 
deviations in the jet energies are added in quadrature. 

The systematic uncertainties in the jet Et scale from the sources 



(a) and (b), as well as from the UE correction are shown in Fig. 5.1 . The 
total systematic uncertainty from these three sources is obtained by adding in 
quadrature the three curves, and is shown as the solid curve in Fig. |5.1| . 

The out-of-cone correction was derived from a Monte Carlo simula- 



tion and accounts for the energy falling outside the jet cone |[49|| . This study 
was done with light quarks; the tt specific corrections take into account dif- 
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ferences with heavy flavor jets. The amount of energy outside the cone of 
R=0.4 is related to emission of low energy gluons from the initial partons, 
and is referred to as "soft gluon" radiation. The correction factor, f oc = 
1+OC(0.4)/Pt(0.4), is a function of the jet Pt corrected for all other effects 
and is given by the equation, 



23.0 (1.0- 0.915 e-o oo^Px) 
Joe = 1-U H • {o.z) 



The correction factor is shown in Fig. |5.2j (a). 

The systematic uncertainty on the jet momentum from the OC cor- 
rection originates from the uncertainty in modeling the radiation of low energy 
gluons in parton showers. To estimate this uncertainty, we use W+l jet data 
and a HERWIG Monte Carlo simulation of W+1 jet events to compare the en- 
ergy contained in an annulus with radii of 0.4 and 1.0 around the jet direction. 
We define a variable F, 

= Px(l.O) - Pt(0-4) 

P t (0.4) V ; 

where Pt(0.4) and Pt(1.0) are the jet momenta corrected using the corrections 
described above (note that Eq. 5.2 is used for R=0.4; for R=1.0 the correction 
is much smaller). The quantity F is the fractional difference of the momentum 
in an annulus with radii between 0.4 and 1.0, calculated for each event using 
the calorimeter towers in that annulus or using the average OC correction. A 
comparison of data and Monte Carlo tests the agreement between the Monte 
Carlo soft gluon radiation modeling and what is observed in the data in that 
annulus. Figure 5~3 shows the mean value of F as a function of the corrected 
Pt (corrected using a cone size of 0.4) for data and Monte Carlo. There is 
a clear difference between the two distributions. This implies that the jet 
shapes in data and Monte Carlo disagree at the few % level. The difference 
between HERWIG and data is shown in Fig. [5.4|. We take this difference as 
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the uncertainty on the out-of-cone correction. Its effect on the top quark 
mass measurement is referred to as the systematic uncertainty from soft gluon 
radiation. 

Similar distributions have been obtained for other sets of data, namely 
Z + 1 jet data and jet data with two 5-tagged jets. Since the statistics for 
the latter sets of data are low, only the H^+l jet data are used. A fit to 
the points of Fig. |5.4| gives a maximum (upper dotted curve) uncertainty of 
(5P T /P T = exp (2.467 - 0.074P T ) + 1.438 (in %). It can be seen that for jets 
typical of those produced in tt events (~ 30 — 90 GeV for M top = 175 GeV/c 2 ), 
the difference between HERWIG and data is < 2%. For softer jets, the differ- 
ence is closer to 4%. 

The systematic uncertainty assigned to the soft gluon radiation ac- 
counts for differences in the energy contained in the annulus 0.4 < R < 1.0 
between data and the Monte Carlo simulation. For the additional energy which 
falls outside a cone of 1.0, we assign an uncertainty of ± 1 GeV. We refer to 
this energy as "splash-out" . 

In summary, Figure |5.2| shows some of the flavor-independent jet 
corrections and their Pt dependence. Figure |5.2| (a) shows the absolute and 
out-of-cone correction factors as a function of the corrected jet Pt- They vary 
from wl.3 at P T = 15 GeV/c to «1.12 for P T > 100 GeV/c. Figure 0(b) 
shows the ratio of the fully corrected jet Pt (P™ 1 ) to the raw jet Pt (P T aw ) 
as a function of the fully corrected jet Pt- Jets from ti events typically have 
a Pt of ~ 30 — 90 GeV/c, for which the average jet correction factor is ~ 
1.45. Figure |5.2| (c) shows the correction factor as a function of P™ w . Finally, 
Fig. |5.2j (d) shows the fraction of momentum measured in the detector before 



the jet corrections as a function of the corrected jet Pt- Figure |5.5| shows the 
overall systematic uncertainty as a function of the corrected Pt of the jets. 
In the 30-90 GeV/c range, the systematic uncertainty on jet energies is about 
4%. 
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Figure 5.3: Fractional difference in corrected jet Pt obtained using cone radii 
of 0.4 and 1.0 as a function of the corrected jet Pt from W+l jet events. The 
circles are the results from the data sample and the triangles are from a sample 
of HERWIG Monte Carlo events which have been processed through the CDF 
detector simulation. 
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Figure 5.4: Systematic uncertainty on the out-of-cone correction as obtained 
in jet events. The abscissa is the fully corrected jet Pt using a clustering 
cone of 0.4. The vertical axis is the difference between data and Monte Carlo of 
the variable F described in the text. The full curve represents a fit through the 
data points, the dotted curves were obtained using the one standard deviation 
values of the fit parameters. 
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Figure 5.5: Total systematic uncertainty on jets as a function of the corrected 
jet P T . 
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5.1.2 Checks on the jet E x scale and its uncertainty 

The procedures used to obtain the jet corrections and their system- 
atics have been checked by applying them to Z + 1 jet events, where the 
Z-boson decays into either e + e~ or fi + fi~. The energy scale for electrons and 
the momentum scale for muons are known to a precision of 0.14% and 0.065% 
respectively In the absence of initial state radiation, such events are 

expected to have zero net transverse momentum. The jet in each event is 
corrected according to the previous prescription, and the quantity 

P T (Z)-P T (jet) 

Fb = p^z) (5 ' 4) 

is calculated, where P^(Z) is in the range 30-150 GeV/c. The lower limit was 
chosen to avoid biases due to the sample selection. The jet recoiling against 
the Z boson is required to have an uncorrected Ex > 8 GeV and \rj\ < 2.4. 
To test the jet energy scale we need a clean environment, i.e., events in which 
there is only one jet recoiling against the Z boson. We therefore require that 
any additional calorimeter cluster have an uncorrected energy Et < 6.0 GeV 
(at any rj). 

To separate detector effects from those due to gluon radiation in the 



initial state, we use the component analysis first suggested in Ref |pT| . We 
compute the direction of the bisector between the Z and the jet directions 
in the transverse plane. The "parallel component" of Ff, is then defined to 
be the component perpendicular to the bisector. Balancing the jet against 
the Z along this component will give information about the jet energy scale. 
Figure |5.6| shows the distribution of this component of Fj, in Z + 1 jet events for 
data and Monte Carlo. The difference in the medians of the two distributions 
is: 

(AF 6 )|| = (3.2 ± 1.5(stat) ± A.l(syst))%. (5.5) 

The 4.1% systematic uncertainty was calculated using the jet energy uncertain- 
ties discussed in the previous section. We conclude that any possible energy 
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scale shift detected by this check is compatible with zero within the evaluated 
uncertainties. 

5.1.3 Jet momentum corrections for tt events 

The tt specific jet momentum corrections are designed to make an 
average correction to the jet momenta to obtain an estimate of the original 
parton momenta 0. The Pt spectra of partons from HERWIG generated tt 



events which pass our experimental selection cuts are shown in Fig. |5.7| . The tt- 
specific corrections account for (a) the difference in the Pt spectrum between 
top induced jets and the flat spectrum used to derive the flavor-independent 
corrections, (b) the energy lost through semileptonic b and c-hadron decays, 
and (c) the multi-jet final state of tt events as compared to di-jet final state 
used to derive the flavor-independent corrections. The correction for these 
three effects are derived using the HERWIG Monte Carlo generator with an 
input top quark mass of 170 GeV/c 2 . The generated events are processed 
using the CDF simulation and reconstructed in the same way as the data 
sample. An average correction factor is determined by first matching (in rj — <p 
space) the reconstructed jets with the generated partons, and then comparing 
the reconstructed jet Pt (after the flavor-independent corrections) with the 
original parton P?. The correction is given by the median of the distribution 
of A = (P^iparton) — P^(jet)) /P^(jet) . This is done as a function of the 
reconstructed jet Pt- 

Figure |5.8| shows the size of the tt-specific correction factors for four 
types of jets: (A) jets from hadronic W decays, (B) average b jets ( no selection 
on decay mode) (C) b jets containing an electron, and (D) b jets containing a 
muon. The general shape of each curve is primarily a result of the difference 
between using a flat jet Pt spectrum and the spectrum appropriate for top 
decays. In particular, this difference is responsible for the rising values of 
the curves at low Pt, and the asymptotic values at large Pt- The larger 



48 



60 - 



C/5 



£40 
> 



20 - 



_ Data 

.... Monte Carlo 



Q L - l -i- -i - J - i- -I i ] i i i I i i i I i i i I i i i I i i i I i i i I i Lj^J i i l 

-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 



(AF b )„ = [(P T (Z)-P T (Jet))/P T (Z)] 



TWJ|| 



Figure 5.6: Parallel component of transverse momentum imbalance between 
the Z and the jet in reconstructed Z + 1 jet candidate events. Both data 
(solid) and Monte Carlo (dashed) are shown (see text). 
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corrections applied to the b jets with a soft lepton are a consequence of the 
amount of energy carried off by undetected neutrinos, and, for jets containing 
a b — > \wX decay, of the fact that muons deposit only ^2 GeV, on average, 
in the calorimeter. 

The flavor-independent and tt-specific corrections bring the median 
reconstructed jet Pt into agreement with the initial parton Pt in ti events. 
The uncertainty on the jet Pt after these corrections is given by the a of the 
A distribution, defined as one half of the separation of the 16th and 84th per- 
centiles of the distribution. For each bin of reconstructed jet Pt, we obtain 
the a of the A distribution, which is then parametrized as a function of the 
reconstructed jet Pt- These uncertainties are shown in Fig. |5.9| for jets from W 
decay and b jets. As above, we display curves for generic b jets (no selection on 
decay mode), for jets containing an electron, and for jets containing a muon. 
These jet Pt uncertainties are input into the kinematic mass fitter (see Sec- 
tion ||) and dictate how much the jet energies can be altered to accommodate 
the applied constraints. 

The jet corrections described above are applied only to the four high- 
est Pt jets in the event, which are assumed to be daughters of the t and t 
decays. Any additional jets beyond the leading four jets are corrected only 
with the "flavor-independent" corrections (excluding the out-of-cone correc- 
tions, see Section |5\2| ) and are assigned an uncertainty of O.IPt © 1 GeV/c. 



This curve is also shown in Fig. 5.9. 



5.1.4 Summary of systematic uncertainties on jet en- 
ergy measurements 

A number of corrections are performed to estimate the original parton 
momenta from the observed jets. The jet energy scale uncertainty is evaluated 
from the uncertainties in the corrections for calorimeter stability, multiple 
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Figure 5.7: Pt distributions for partons from top quark decays obtained from 
the HERWIG Monte Carlo program after simulation of detector response and 
including the effects of the Top Mass Candidate Sample data selection. The 
solid line indicates the distribution for light quarks from the W — > qq' decay 
and the dashed line is the distribution for b quarks. 
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Figure 5.8: The tt-specific corrections applied to jets according to available b- 
jet information. The curves show the fractional change to the corrected jet Pt 
after all "flavor-independent" jet corrections have been applied. The curves 
are for: (A) jets from the decay of W bosons, (B) jets from all b quarks (no 
selection on decay mode), (C) jets from b quarks containing an electron, and 
(D) jets from b quarks containing a muon. 
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Jet P T (GeV/c) 



Figure 5.9: Fractional uncertainty in the estimated parton Pt as a function 
of the jet Pt after the flavor-independent jet corrections are applied. The 
uncertainty shown on the vertical axis is given as a fraction of the jet Pt- 
Curves (A) through (D) have the same meaning as for the previous figure. 
The curve labeled (E) is used for the jets beyond the four highest-Px jets and 
is applied only to the Pt within the cone of radius R=0.4. 
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interactions, calorimeter response, jet fragmentation, underlying event, out-of- 
cone correction, and splash out. Figure [575] shows the dependence of the overall 
jet energy scale uncertainty on the corrected jet Pt- The total systematic 
uncertainty varies between 7% for jets with corrected Pt of 20 GeV/c and 
3.5% for jets with P T =150 GeV/c. 

We do not assign a separate systematic uncertainty to the top specific 
corrections. Such uncertainties may arise from modeling of initial and final 
state gluon radiation, and modeling of the primary parton collision. We discuss 
these uncertainties in Section || 



5.2 Measurement of other calorimeter variables 

To measure the top mass we apply energy-momentum conservation 
to the process pp — ► tt + X, with subsequent decay of the t (t) into W + b (b) 
(see Section |||). Here, X is the unspecified particles which recoil against the 
tt system. The calorimeter provides the measurement of X?, the transverse 
momentum of X. The quantity Xt is computed from the energy left over 
after the lepton and the four jets from the tt system are removed from the 
total measured energy. This leaves two terms: 

Njets 

X T = U T +Y. Mjet) (5.6) 



i=5 



Each component of the unclustered energy, Ut, is defined as the vector sum 
of the energies in the calorimeter towers after excluding the primary lepton 
and all the jets with raw Et > 8 GeV and \r]\ < 3.4 in the event. Using a tt 
Monte Carlo (M top =175 GeV/c 2 ) we find a distribution in U x with (U x ) ~ 
and a = 15.8 GeV for events which enter into the mass subsamples. The same 
distribution for the data has a mean consistent with and a a = 14.9 GeV. 
The Monte Carlo and data distributions in U x are shown in Fig. |5.10| . Similar 
results are obtained for the y component. 
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Figure 5.10: Unclustered energy for the mass sample used here and for tt 
Monte Carlo (M top =175 GeV/c 2 ). Only the x component is shown. 
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Each component of the unclustered energy is corrected with a single 
factor /u. e . = 1.6, based on studies of the recoil energy that the calorimeter 
measures in Z-boson events with no extra jets, where the Z boson is well 
measured by the two leptons it decays into fl52| . The final mass value is not 
sensitive to the value of / u , e .. For example, using f u e ^ = 2.0 makes a negligible 
change in the reconstructed top quark mass (0.2 GeV/c 2 ), hence it is not 
included in the Table of systematics in Section |9.7| . The uncertainty with which 
each component of Ut is measured is taken to be 100% (added in quadrature 
to 1 GeV). The jets beyond the four with the highest Et are corrected only 
within the cone of 0.4, so as to avoid counting the out of cone energy twice (it 
is already included in the unclustered energy). The uncertainties on these jet 
energies were discussed in Section [5.1.3| and shown in Figure 5J3. 



Another quantity that can be estimated from the calorimeter mea- 
surement is the E T . It is calculated using the following expression: 

_ 4 

-% T = E T (lepton) + E T (jet) + X T (5.7) 

The above expression shows that the Ex measurement is highly correlated 
with the jet energy measurements, and therefore it is not considered as an 
independent measurement in the mass fitting. As discussed in Section |], it is 
only used as a starting value for the neutrino's transverse momentum when 
the overall mass fit is performed. 



Section 6 
MASS FITTING 



The kinematics of events in the decay channel pp — > tt — > iuqq'bbX 
are over-constrained by the number of measured quantities and the number of 
applicable energy-momentum conservation equations of production and decay. 
This allows for complete reconstruction of the four-momenta of the particles 
in the decay chain and hence an event-by-event top mass determination. 

In this section we discuss the methods used for event reconstruction 
and then study the validity of the algorithms using tt Monte Carlo events. 
Effects due to combinatorics, wrong parton assignments and shapes of back- 
grounds on the top mass measurement are also discussed. 



6.1 Event reconstruction 

The first step in the reconstruction is the estimation of the four- 
momenta of the decay products of the tt pair: the lepton and the four jets. 
Electron and muon measurements, resolutions and identification are discussed 



in section |3J]. The four leading jets, as defined in Section are assumed to 
be the q, q, b, and b quarks from the tt decay chain. According to a HERWIG 
Monte Carlo plus detector simulation, this assumption is correct 55-72% of the 
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time, depending on the number and type (SVX or SLT) of tags (see Table |Q| ). 
The momenta of the reconstructed jets are corrected as described in Section [5[ 
The direction (i.e., r\ and 4>) of each parton is assumed to be the same as the 
direction of the associated jet. The masses of the partons are assumed to be 
0.5 GeV/c 2 , except for h and b quarks which are assigned a mass of 5.0 GeV/c 2 . 
The resolutions on the jet energy measurements are discussed in Section [5[ 

The mass fitting algorithm applies the constraints implied by the 
production and decay of a ti pair to evaluate an event-by-event mass. The 
hypothesis of standard model ti implies the production process 

pp-^t + i+X, (6.1) 

followed by the decays 



t -> W + + b, (6.2) 

i -> W~ + b, (6.3) 

W ± -»■ £ ± + u, (6.4) 

W T -> q + q'. (6.5) 

The quantity X, in pp — > t+i+X, represents the unspecified particles 
recoiling against the ti system. Only two components of X are measured, as 
discussed in Section fT2[ 

An estimate of the top quark mass is obtained on an event-by-event 
basis after minimizing a x 2 - m general, the \ 2 definition is not unique, in that 
any formulation which expresses the constraints implied by the measurements 
and four-momentum conservation is equally valid. We have chosen a particular 
formulation of the x 2 ; which is minimized using the program MINUIT |5v|. An 



alternate method, the SQUAW kinematic fit [54], has also been used and the 



results are essentially identical. We describe here both of these fitters. 
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6.1.1 Mass fitting using MINUIT 

The x 2 expression which uses the MINUIT minimization routines ap- 
plies energy and momentum constraints to the above production and decay 
chain to obtain six effective constraints: (1,2) the two transverse momentum 
components of the ti + X system must be zero, (3) the invariant mass of the 
Iv system must equal the IT-boson mass, M\y, (4) the invariant mass of the 
qq' system must equal Mw, and (5,6) the two three-body invariant masses 
must each equal the top quark mass, M t . The relevant unmeasured quantities 
are then the three momentum components of the neutrino and the top quark 
mass. The system may therefore be solved by minimizing a two-constraint 
chisquare. The chisquare expression used to obtain the present results is: 

2 / ~ ,\ 2 



2 = v (Pt-Pt) (Uj-Ul) {M tv -M w f 

X. / i 2 / j 2 2 

£,jets a P T i=x,y IT U' i ' M w 

| (M n -M w ) 2 | {M lv3 -M t f | (M^-M f f 

°M W a M t a M t 

The notation is as follows: £ signifies the primary lepton in the event, v refers 
to the inferred neutrino, and j refers to one of the four leading jets in the event. 
The first sum is over the primary lepton and all jets with raw Et > 8 GeV 
and \tj\ < 2.4. The second sum is over the transverse components of the 



unclustered energy p5| , discussed in Section [5T2| , plus those of the energies of 
jets with 2.4 < \r)\ < 3.4. The hatted symbols in the sums represent quantities 
altered by the fit procedure, whereas unhatted symbols represent the input 
values. The uncertainties on the energy of the primary lepton, the jets, and 
the unclustered energy are discussed in Sections |34] and [5|. The W-boson 
mass, Mw, is taken to be 80.4 GeV/c 2 f3"Ifl , <tm w is set to 2.1 GeV/c 2 f5"5| , and 
is set to 2.5 GeV/c 2 . The results are insensitive to the values used for o"m w 
and out- The quantity Mt v is the invariant mass of the primary lepton and 
the neutrino and Mg U j is the invariant mass of the primary lepton, neutrino, 
and one of the four leading jets. Of the remaining three jets, we assign two to 
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the decay products of the W boson in order to calculate Mjj. The third jet 
is then combined with the other two jets to form the three body mass Mjjj. 



The issue of combinatorics is discussed in Section |T2 . 

The first two constraints are that the total transverse momentum 
components of the ti + X system are zero. These constraints are imposed by 
setting the neutrino transverse momenta to exactly balance the sum of the 
current Pt and U T values. The other four constraints appear as explicit terms 
in the \ 2 - This x 2 yields two minima which correspond to the two solutions 
for the neutrino longitudinal momentum in the W decay. This is referred to as 
the P£ ambiguity. After minimizing this \ 2 with respect to the collective set of 
transverse momenta, Pt, for the jets and the charged lepton, the unclustered 
energy, Ut, the ^-component of the neutrino momentum, and the top mass, 
M t , for the event, we obtain an event-by-event determination of the top quark 
mass. 



6.1.2 The SQUAW fitter 

The SQUAW fitter is a general kinematic fitting program that can be 
used for any production and decay processes, provided that there are enough 
constraints [|I], It has been used to measure the top quark mass in the 
lepton+jets channel and for the all-hadronic decay channel ^ . 

In brief, it applies energy-momentum conservation to the five pro- 



cesses (|6.1]) - (|6Tq) , thus providing 20 equations, i.e., 20 constraints, for the 
measured quantities and their uncertainties. It uses the measured W mass, 
M\y = 80.4 GeV/c 2 . In the fit an uncertainty is assigned to the W mass in 
order to take into account the expected W width of 2.1 GeV/c 2 . Additional 
ingredients of the kinematic fit: 

• The measured quantities are: the lepton, the four leading jet momenta, 
and Xt- 
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• For each event there are 18 unknowns. These are: energy and P z of X 
(2), 3- momenta of t and t plus the top mass (7), 3-momenta of the W 
bosons (6), and the 3-momenta of the v (3). 

This is then a 5- vertices, 2-constraints fit, 5V-2C in SQUAW's language. Notice 
that the v momentum is considered an unknown quantity. This is because the 
JZl T is highly correlated with the jet momentum measurements. The calculated 
value of JZl T is used as a starting point to help with the convergence of the fit. 
Lagrange multiplier techniques are used to solve the 20 equations. The final 
X 2 has contributions from all 20 equations. 

One of the differences with the MINUIT algorithm is that SQUAW 
works with the 4-vectors, hence it allows the angles of the lepton and jets to 
vary within their uncertainties. The momentum magnitude and angles are 
assumed to be uncorrelated. 

The results of the two methods for a given event are very close. In 
the 76 event data sample the masses obtained with the two methods (using 
the mass from the lowest x 2 solution in each case) differ on the average by 0.1 
GeV/c 2 , and in 70% of the events the absolute value of the mass differerence 
is less than 0.5 GeV/c 2 . 

6.2 Combinatorics 

There is always some ambiguity in how to assign the four leading jets 
to the four relevant partons. If none of the jets is tagged as a b candidate, by 
either the SVX or SLT algorithm, then there are 12 different ways of assigning 
jets to the b and b partons. Combined with the ambiguity, there are 
then 24 combinations, or configurations, per event. If one jet is tagged as 
a b candidate, we require that it is assigned to a b or b parton, and this 
reduces the number of allowed combinations to 12. If two jets are b tagged, 
there are four combinations. Of the above combinations the solution with the 
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lowest x 2 is chosen, and that solution is required to have x 2 < 10. The latter 
requirement defines criterion 9 of the Top Mass Candidate Sample described 
in Section p.3| . We have not found a satisfactory method for improving the 
top mass resolution by including any solutions with x 2 values larger than the 
lowest one, and therefore we take the lowest x 2 solution as the best estimate 
of the top mass for each event. 



6.3 Impact of gluon radiation 

A substantial fraction of ti events are expected to contain extra jets 
resulting from gluon radiation. From a HERWIG Monte Carlo plus detector 
simulation, we find that ~40% of events have one or more jets which do not 
correspond to the partons from the ti decay. These extra jets may be produced 
during the production of the tt pair (initial state radiation) or in the decay 



stage (final state radiation) |56|. From a theoretical perspective, whether 
or not the extra jet(s) are to be included in the fit depends on whether the 
gluon was radiated during production of the ti pair or during its decay. If the 
radiation comes from the production stage, then it should not be included in 
the mass fit. If the radiation is produced from a quark in the decay stage, 
then it should be included as one of the decay products [ 155) . 

From an experimental perspective, the radiation results in jets which 
may or may not have been produced in the ti decay process. On an event- 
by-event basis, production and decay stage radiation cannot be differentiated 
from each other or, for that matter, from the partons from the ti decay (unless 
the jet is b tagged). Gluon jets which come from decay stage radiation are more 
correlated with the partons emerging from the hard scatter, and therefore one 
can consider merging jets which are close in i]—<j) space. It is also possible to try 
all unique permutations of four jets among all the reconstructed jets. However, 
taking a fifth jet into consideration increases the number of combinations by a 
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factor of 3, 4, and 5 for the 2, 1, and fe-tag cases respectively. This increase 
in the number of solutions reduces the probability for choosing the correct 
jet assignment. The mass reconstruction presented here does not implement 
either of these possibilities. Our approach is to assume the model of initial and 
final state radiation in the Monte Carlo simulation is correct, and to associate 
a systematic uncertainty with this assumption. 

6.4 Results of the kinematic fit on simulated 

ti events 

The reconstructed-mass distribution obtained by fitting simulated ti 
events depends on the intrinsic resolution of the detector, and, more impor- 
tantly, the ability to correctly associate the daughter partons from a ti decay 
with the observed jets. Both combinatorics and gluon radiation play a role in 
degrading the resolution of the top quark mass measurement. In this section, 
we discuss the performance of the mass fitter by dividing events (which enter 
into one of the four mass subsamples) into three categories: 

1. Correctly Assigned Events: Each of the four leading jets are within 
AR < 0.4 of a parton from the ti decay and are correctly associated with 
the appropriate quark by the lowest y 2 solution satisfying any imposed 
tagging requirements. The jet-parton match is required to be unique. 

2. Incorrectly Assigned Events: Each of the four leading jets are within 
AR < 0.4 of a parton from the ti decay and each jet-parton match is 
unique, but the configuration with the lowest x 2 is n °t the correct one. 

3. Ill-Defined Events: The four leading partons from the ti decay cannot 
be uniquely matched (AR < 0.4) to the four leading jets in the event. 
Such events often have extra jets produced from either initial state or 
final state radiation. 
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The fractions of events falling into each of these categories are estimated using 
a HERWIG tt Monte Carlo plus detector simulation. These fractions depend 
on the 6-tagging information in the event. For example, having two 6-tagged 
jets in an event reduces the probability that one (or more) of the leading 
four jets is a gluon jet. The fractions of events falling into categories 1-3 
above, and the width of the reconstructed-mass distribution for each of the 
four mass subsamples are shown in Table |6.1j . The widths are calculated as 
half the difference between the 16 th and 84 th percentiles of the reconstructed- 
mass distributions. The reconstructed-mass distributions for the four mass 



subsamples are shown in Fig. 5.1 





Event fractions (%) 


Width 


Data sample 


1 


2 


3 


(GeV/c 2 ) 


SVX Double 


49 ±2 


23 ±2 


28 ±2 


19.9 


SVX Single 


30 ±1 


26 ±1 


44 ± 1 


24.2 


SLT 


26 ±2 


31 ±2 


43 ±2 


25.0 


No Tags 


23 ±1 


32 ±1 


45 ±1 


26.9 



Table 6.1: Fractions of tt events falling into categories 1-3 described in the text. 
The last column shows the width of the distribution of reconstructed masses 
for each subsample. The width is taken to be half the difference between the 
16 th and 84 th percentiles of the relevant mass distribution. 



As Table |TT] shows, the fraction of correctly assigned jets increases as 
the number and purity of b tags increase. Figure |T| shows the reconstructed- 
mass distributions for events in each of these three categories. When the 
correct jet-parton assignments are made (category 1), the resolution is ~ 
13 GeV/c 2 , while for categories 2 and 3 it is ~36 and 34 GeV/c 2 respectively. 
As Fig. |6.2| demonstrates, the mass resolution is dominated by incorrect as- 
signment of jets to partons from the tt decay. For Double SVX tagged events, 
where nearly half of the events have the four leading jets correctly assigned 
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to the tt decay products, we obtain the best resolution on the reconstructed 
mass. 

A priori, it is not obvious whether events which have the jets mis- 
assigned to the tt daughter partons contain information on the top quark 
mass. This is quantified by studying the sensitivity of the distribution of 
reconstructed masses to changes in the input value of the top quark mass. 
We examine the events in categories 1-3 separately in order to determine if 
the misassigned events contribute information to the top quark mass mea- 
surement. For each category of events, we evaluate the rate of change of the 
median of the reconstructed-mass distribution as we vary the input value of 
the top quark mass. Larger changes in the median imply greater sensitivity 
to the top quark mass. Figure shows the median reconstructed-mass as a 
function of the input top quark mass. Events from all four mass subsamples 
are included in the distributions. The four distributions correspond to events 
in (a) category 1, (b) category 2, (c) category 3, and (d) the three categories 
combined. We find that the events in which the jets are correctly assigned 
to the partons have the largest slope (0.90), while incorrectly assigned events 
have a slope of 0.62 and ill-defined events have a slope of 0.48. Correctly as- 
signed events (category 1) do not have a slope of 1.0 because the top-specific 
corrections (see Section |5.1.3j ) are derived using a specific input top quark 
mass of 170 GeV/c 2 . We conclude that the events with incorrect jet-to-parton 
assignments do in fact contain information on the top quark mass, since the 
slope is not zero. However, because of the smaller slope and larger width of 
the reconstructed-mass distribution, incorrect combinations degrade the res- 
olution of the top quark mass measurement. The slopes for each of the four 
subsamples in each category are shown in Table |6.2| . The slopes vary from a 
maximum of 0.81 for SVX Double tags to a minimum of 0.62 for SLT tagged 
events. Since SVX double-tagged events have the largest slope, narrowest 
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width and lowest background, they generally yield the best precision on the 
top quark mass measurement (for equal size subsamples). 



Slopes 

Data sample 12 3 Combined 



SVX Double 
SVX Single 
SLT 

No Tags 



0.89 ±0.02 
0.90 ±0.01 
0.93 ±0.02 
0.90 ±0.01 



0.48 ±0.08 
0.60 ±0.04 
0.68 ±0.05 
0.62 ±0.03 



0.57 ±0.05 
0.52 ±0.02 
0.38 ±0.04 
0.47 ±0.03 



0.81 ±0.03 
0.72 ±0.02 
0.62 ±0.03 
0.68 ±0.02 



Table 6.2: Rate of change ("slope") of the median reconstructed mass with 
the input value of the top quark mass, for the four mass subsamples. For each 
subsample, we show the slope for the three categories of events defined in the 
text, both separately and combined. 



6.4.1 Mass reconstruction in other tt decay channels 

Although the fitting procedure assumes that the candidate tt events 
have decayed through the tt — > (e or fi)uqq'bb channel, there is a non- negligible 
contribution from top events decaying through other channels. The additional 
acceptance from other decay channels comes mostly from events where either 
an electron or a r from the t or i decay is misconstrued as a jet or from events 
with a leptonically decaying r. In either case, two b jets are still present. 
The fourth jet can be produced through gluon radiation. Table O gives the 
expected contributions of various decay channels to the candidate tt sample 
and to the subsample with at least one SVX or SLT tag. It shows an 11% 
contribution from r events and 4% contribution from ee, e/x and fifi events. 



Figure I^J shows the reconstructed-mass distribution for events from 
these decay channels. The inset shows how the median of the reconstructed- 
mass distribution changes with the input value of the top quark mass used in 
the simulation. The relatively low, but non-zero value of the slope indicates 
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that these events also provide information about the top quark mass. The 
signal templates, to be discussed in Section [7|, include contributions from these 
channels, hence we do not expect any bias on the fitted top mass from these 
events. 



Channel 


Top Mass Candidate Sample 


Tagged Events 


e+jets 


0.423 ±0.008 


0.424 ±0.010 


/x+jets 


0.426 ±0.008 


0.430 ±0.010 


e — r(had) 


0.017 ±0.002 


0.017 ±0.003 


e — r(lep) 


0.007 ±0.001 


0.007 ±0.002 


ee 


0.012 ±0.002 


0.011 ±0.002 


e — ii 


0.023 ±0.002 


0.024 ±0.003 


\l — r(had) 


0.017 ±0.002 


0.017 ±0.003 


\i — r(lep) 


0.005 ±0.001 


0.006 ±0.002 


n — n 


0.004 ±0.001 


0.004 ±0.001 


r(had) + rilep) 


0.002 ±0.001 




r(/ep)+jets 


0.063 ± 0.004 


0.058 ±0.005 



Table 6.3: Fractional contribution (according to Monte Carlo simulation) of 
lepton±jets events to tt events in the Top Mass Candidate Sample and the 
Tagged subsample. The Tagged subsample includes events with at least one 
6-tagged jet. Similar numbers are found for other subsamples {i.e. SVX Dou- 
ble, SVX Single and SLT). lep and had denote leptonic and hadronic decays, 
respectively, for the r lepton. A dash indicates that no events were found in 
the category. 



6.5 Mass reconstruction in won-tt events 

Non-tt events are also present in the data samples. For all the samples 
considered, the dominant background is expected to be from production of W 
bosons in association with extra jets. The background shape is modeled with 
the VECBOS Monte Carlo simulation. As with tt events, we fit the background 
events using the x 2 defined in Section |6|. Since the sample of events does not 
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contain tt, one does not expect any resonant peaks in the reconstructed-mass 
spectra. The reconstructed-mass spectrum for VECBOS events which have at 



least one 6-tagged jet is shown in Fig. This distribution is compared to 
the distributions for tt events with input top quark masses of 140, 175, and 
200 GeV/c 2 . It is observed that for a top quark mass of 140 GeV/c 2 , the 
signal and background peak at nearly the same value of reconstructed mass. 
However, the ti events are more sharply peaked than background, and therefore 
there is still shape discrimination between the two. As the top quark mass 
increases, the reconstructed-mass distribution for tt events is clearly separated 
from the background. Since we include a background constraint in the top 
quark mass likelihood fit (see Section [7|), differences in shape between signal 
and background events are not required. However, the shape differences do 
improve the resolution on the top quark mass measurement. 



6.5.1 Comparisons of VECBOS with data 

The background modeling is checked by comparing the reconstructed 
W+A jet mass distributions from some data samples with the appropriate 
distributions from the VECBOS simulation. The data samples consist of events 
that fail only one of the top sample criteria. The samples compared are: 

• events failing only the lepton isolation criteria; 

• events having fewer than three jets with Ex > 15 GeV and \r)\ < 2 (see 
Section |3^3| ); 

• events with a non-central (1.1 < \rj\ < 2.4) primary electron. 



Each of these samples fails one and only one of the top sample criteria. Fig- 
ure |6.6| shows the reconstructed-mass spectrum for candidate W+4 jet events 
in which the primary electron is not isolated from jet activity in the event. The 
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requirement that the lepton is non-isolated makes it more likely that the se- 
lected data events are from multi-jet or bb production. The data are compared 
to the distribution from VECBOS events which also failed the lepton isolation 
criteria. The similarity of the two distributions shows that the non-W/Z com- 
ponent of the background is well modeled by the VECBOS simulation. The 
fraction of tt events in the data sample is expected to be ~9%. A Kolmogorov- 
Smirnov test applied to these two distributions yields a 36% confidence level 
for agreement. Figure |6.7| shows a similar comparison in which the events 
are required to have no more than two jets with Et greater than 15 GeV. 
This sample has an estimated tt contribution of about 0.7%. A Kolmogorov- 
Smirnov test applied to these two distributions yields a 45% confidence level 



for agreement. Figure |6J| compares samples of events in which the primary 
electron was reconstructed in the PEM (1.1 < \rj\ < 2.4). We expect little or 
no dependence of the reconstructed mass on the r\ value of the primary elec- 
tron, as evidenced by the similarity between this VECBOS distribution and 
the one in Fig. |6.5| . This sample is estimated to have a tt fraction of 0.2%. 
A Kolmogorov-Smirnov test applied to these two distributions yields a 33% 
confidence level for agreement. 

We expect the events in these three data samples to be predominantly 
from the same sources as described in Section [4.3| , but in different proportions. 
In all three cases the VECBOS simulation agrees with the reconstructed-mass 
distribution in the data. Therefore we assume that the VECBOS simulation 
models satisfactorily the reconstructed-mass distribution of the background 
events in the mass subsamples. 
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Figure 6.1: Reconstructed-mass distributions for HERWIG tt events 
(M top =175 GeV/c 2 ) for the mass subsamples: (a) SVX Double, (b) SVX 
Single, (c) SLT (no SVX), and (d) No Tags. The 16 th and 84 th percentiles for 
each distribution are indicated by the arrows on the figures along with their 
values. 
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Figure 6.2: Reconstructed mass for M top = 175 GeV/c 2 tt events which enter 
into any of the four subsamples. The black filled histogram shows the distri- 
bution for those events for which the selected jet-parton configuration was also 
the correct one (category 1). The lightly shaded histogram shows the distri- 
butions for which a correct assignment could be defined, but was not selected 
(category 2). The darker shaded histogram shows the distribution for events 
where a correct assignment was ill-defined (category 3). The solid line shows 
the three distributions combined. 
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Figure 6.3: Median of the reconstructed-mass distribution as a function of the 
input top quark mass used in the simulation. The medians are evaluated from 
mass distributions which include events from all four mass subsamples. The 
figure demonstrates the relative sensitivity of the reconstructed-mass distri- 
bution to the input top quark mass. We show the distributions for events in 
(a) category 1, (b) category 2, (c) category 3, and (d) the three categories 
combined. The slopes indicated were evaluated using a linear fit (dashed line) 
to the data points. 
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Figure 6.4: Reconstructed masses for tt events (M top =175 GeV/c 2 ) which 
decay into decay channels other than e+jets or /x+jets, but fit the hypothesized 
tt — > (e or ji)uqq'bb decay chain (points). Most of these events are due to W 
decays into r leptons (see Table 6.3). Events from all four mass subsamples 
are included. For comparison, the dashed distribution corresponds to events 
decaying through the tt — > (e or fi)uqq'bb channel, but the lowest x 2 solution is 
incorrect (i.e., events in categories 2 and 3). The inset shows how the median 
of the reconstructed-mass distribution changes with the input value of the top 
quark mass used in the simulation. 
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Figure 6.5: Reconstructed-mass distribution for background events from the 
VECBOS Monte Carlo simulation (points). Also shown are the reconstructed- 
mass distributions from the HERWIG Monte Carlo simulation using input top 
quark masses of 140, 175, and 200 GeV/c 2 . In all cases, events are required to 
have at least one SVX or SLT tagged jet. Each distribution is normalized to 
have unit area. 
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Figure 6.6: Reconstructed-mass distribution for events containing at least four 
jets, FJ T , and a non-isolated lepton. The expected fraction of tt in this sample 
is ~9%. The points are data and the histogram is the VECBOS distribution. 
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Figure 6.7: Reconstructed-mass distributions for events with an isolated lep- 
ton, JZl T , and no more than two jets with Et > 15 GeV. The expected fraction 
of ti in this sample is ~0.7%. The points are data and the histogram is the 
distribution from VECBOS. 
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Figure 6.8: Reconstructed-mass distribution for events containing at least four 
jets, F^t , and a non-central electron (reconstructed in the region 1.1 < \r)\ < 
2.4). The expected fraction of ti in this sample is ~0.2%. The points are data 
and the histogram is the VECBOS distribution. 
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DESCRIPTION OF THE 
LIKELIHOOD PROCEDURE 

A likelihood procedure is used to extract a measurement of the top 
quark mass from the reconstructed-mass distributions of the data samples 
and the tt signal and background models, along with the constraint on the 
background fractions. This section describes the likelihood and discusses its 
validation with the help of simulated experiments. 

7.1 Parametrization of the reconstructed-mass 
distributions 

We use the HERWIG Monte Carlo plus detector simulation to model 
the shape of the reconstructed-mass distribution in tt events. Event samples 
are generated at several different values of the top quark mass ranging from 
120 to 220 GeV/c 2 . The VECBOS Monte Carlo program is used to model 
the shape of the background distribution. The tt and background samples are 
processed using the CDF simulation, and the same analysis is applied to them 
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as to the data sample. Histograms of reconstructed masses from these samples 
are referred to as templates. 

Since we generated templates for only a finite number of input top 
quark masses, extraction of a measured top quark mass from the data will 
require an interpolation. This can be achieved in one of two ways. The first 
consists of fitting the data separately at each input top quark mass value to 
a combination of signal (at the given mass) and background. The maximum 
likelihood is then interpolated from the resulting likelihood values. The second 
method requires that the templates themselves be interpolated as a function 
of input top quark mass. The signal templates then become a function of 
both the input top quark mass and the reconstructed mass. The background 
templates are only functions of the reconstructed mass. The likelihood can 
then be defined as a smooth function of both input mass and reconstructed 
mass, and no further interpolation is needed. Previous publications 0, |J, 
with lower statistics, have used the first method of interpolation. However, 
the second method employs a single interpolation process and uses optimally 
the finite Monte Carlo statistics in the templates. We have adopted the latter 
method for this analysis. 

7.1.1 Signal parametrization 

A single function, f s , is used to model the distribution of recon- 
structed top masses for ti events for any given value of the input top quark 
mass between 120 and 220 GeV/c 2 : 

f s (M t , P k ) = N [P 6 fi(M t , P 1)2 , 3 ) + (1 - P 6 ) f 2 (M t , P 4 , 5 )] , (7.1) 

where: 

pi+P 2 

/i(M t ,P 1)2 , 3 ) = r(1 3 + p 2) (M t -P 1 ) ft e- i W-*>, (7.2) 

1 i / M t -P 4 \ 2 

/ 2 (M t ,P 4>5 ) = -^==-e~n * ) , (7.3) 
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Pi = c*i + a i+6 M top . (7.4) 

For each mass subsample, six pairs of parameters (ctj, a i+6 ) are needed to de- 
scribe how the distribution of reconstructed mass (M t ) evolves with the input 
top quark mass (M top ). For example, a 4 and a w (a 5 and an) describe how 
the mean (width) of the Gaussian portion of the reconstructed-mass distribu- 
tion changes with M top . The parameter values and their covariance matrix 
are obtained by a chisquare fit to the templates |^7|. Six of the 18 templates 
for the SVX Single sample are shown in Fig |7.1| together with the predictions 
obtained from the fit parameter values. Figs [7.2| - [7.4| show the same six tem- 
plates for the SVX Double, SLT, and No Tag subsamples respectively. The fit 
chisquares per degree of freedom (dof) are: 1.17 for 555 dof, 1.07 for 335 dof, 
0.96 for 454 dof, and 1.36 for 589 dof, for these four subsamples, respectively. 



7.1.2 Background parametrization 

The fitting of the distribution of reconstructed masses from VEC- 
BOS is performed in a similar fashion to the signal templates, but with fewer 
parameters and no dependence on M top . For the tagged subsamples, the 
background distribution shape can be described by fx, whereas the No Tag 
subsample requires the additional freedom of / 2 to adequately describe its 



shape. Figure |7.5| shows the parametrizations of the background distributions 
for the SVX tagged, SLT (no SVX) tagged, and No Tag events. Because of 
limited statistics and low probability for obtaining two SVX tagged jets in the 
VECBOS Monte Carlo simulation, we assume the same background shape for 
SVX Double and SVX Single tag events. The mass measurement is insensi- 
tive to this assumption because the expected background fraction for double 
tag events is only 4%. In Section |6.5.1| we compared distributions from top- 
depleted data samples with analogous VECBOS distributions to show that the 
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VECBOS Monte Carlo simulation models the shape of the tt backgrounds quite 
well. 



7.2 Definition of the mass likelihood 

The value of M top for each mass subsample is obtained from a max- 
imum likelihood procedure applied to the observed mass distribution. The 
procedure allows the template fit parameters and the background fraction Xb 
to vary about their central values within their respective uncertainties. The 
only parameter which is entirely unconstrained in the mass likelihood fit is 
M top . The reconstructed masses for the events in each of the four mass sub- 
samples are tabulated in Appendix A. Since the subsamples are statistically 
independent, the probability of observing these four sets of masses can be 
expressed as a product of four individual likelihood functions, one for each 
subsample. These four likelihoods have the same form: 



where: 



f-'shape X ^-background ^parami '-5) 



-^events 

C-shape = II [(l-Xb)fs(M t ,M top ,a)+x b f b (M i J)\, (7.6) 



^-background P\p^b)i V " V 

C-param = expj-^ (« - a ) T U~ l (<J - cTq) 

+ 0-fay V - x 0-&j\\. (7.8) 



The likelihood C s ha P e is the joint probability density for a sample of N events 
reconstructed masses Mi to come from a parent distribution with background 
fraction x& and signal fraction (1 — x ). The background likelihood ^-background, 
discussed in Section |4.3.3| , constrains the fraction of background events to the 



expected value within its uncertainties (see Fig. |4.1| ). The expected back- 



ground fraction and number of background events are related via Nb = xy, x 
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N obs , where N obs is the number of observed events for that subsample. The 

— # 

vectors a and (3 determine the shapes of the signal (f s ) and background (fa) 
distributions. They are constrained by C param to agree with the nominal val- 
ues, do and (3q, via their covariance matrices U and V respectively. The 
inclusion of C param in the likelihood definition is due to the finite statistics 
of the Monte Carlo samples used to determine f s and fj,. Furthermore, by 
parametrizing the signal probability f s as a continuous function of M top , the 
likelihood is inherently a continuous function of M top as well. 

To extract the top quark mass for each subsample, we minimize 

— * 

— log£ with respect to M top , x b , a and (3. The statistical uncertainty on 
M top is taken as the change in M top which results in a 0.5 unit increase in 

— log C along the line on which — log C is minimized with respect to variations 
in all the other fit parameters. The statistical uncertainty has contributions 
not only from the finite statistics in the data sample, but also from the un- 
certainty in the expected background and the finite statistics in the mass 
templates. However, the latter two contributions account for less than 1% of 
the total statistical uncertainty. The top quark mass and its statistical un- 
certainty for the four subsamples combined are extracted in the same way as 
above from the product of the four subsample likelihoods. 

7.3 Tests of the likelihood procedure on sim- 
ulated experiments 

The performance of the likelihood scheme was tested using simulated 
events from Monte Carlo programs. We performed a large number of simulated 
experiments, each consisting of four subsamples with the same numbers of 
events (N l ohs) i = 1,...,4) as observed in the four data subsamples. Each 
experiment subsample contained iV£ background events and N l s = (N* bs — N^) 
ti events, where iV£ is a binomial fluctuation of the expected background. 
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The Y<i Nl and Yh distinct mass values for each simulated experiment were 
chosen at random from the discrete templates for signal and background events 
(Figs |7. 1| - |7. 5| ) . The four sets of masses were fit using the same likelihood 
procedure that was used to fit the data sample. Each simulated experiment 
yielded a fitted top quark mass, a statistical uncertainty and a maximum 
likelihood value. The self-consistency of the likelihood procedure was tested 
by comparing these returned values with expectations. 

Figure |7.6j shows the distribution of returned masses from the like- 
lihood fit for input top quark masses of 150, 175, and 200 GeV/c 2 . The 
curves are fits to Gaussians, and have central values of 149.8, 174.8, and 200.2 
GeV/c 2 , and o of 5.8, 6.8, and 7.6 GeV/c 2 . In each case the mean of the 
distribution is consistent with the input value, which demonstrates that the 
procedure introduces little or no bias into the top quark mass measurement. 
The a of the distributions reflects the expected statistical uncertainty on the 
top quark mass measurement for experiments which have the same expected 
background and 6-tag composition as our Run 1 data sample. Based on the 
fitted a's one expects to achieve a statistical uncertainty on M top of ~4%. 

The statistical uncertainty returned by the likelihood procedure should 
reflect the deviation of the returned top quark mass from the input value. The 
pull, defined by 

Pull = Mexp ~™ input , (7.9) 

is used to check the consistency between the measured deviation on the top 
quark mass and the estimated statistical uncertainty. In the above expression, 
M exp is the fitted top mass value returned by the likelihood, M input is the 
input value used to generate the (simulated) experiment, and a^ 1 is the sta- 
tistical uncertainty on M exp returned by the fitter. Figure [7/7| shows the pull 
distribution for the simulated experiments generated for M top = 175 GeV/c 2 . 
The width is close to unity, which indicates that the statistical uncertainty 
returned by the fitter accurately reflects the deviation of the fitted value from 
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the input value. Alternately, in Figure |778| , we take slices in a s ^ at , and evaluate 
the width of the corresponding (M exp — M input ) distribution. The points have 
a slope of 0.92±0.09, which supports using the statistical uncertainty returned 
by the fitter as a measure of the statistical uncertainty for a given experiment. 
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Figure 7.1: Reconstructed- mass distributions from simulated it events for sev- 
eral input values for the top quark mass used in the simulation. The overlayed 
curves are predictions from the parametrization of templates at 18 different 
top mass values. The distributions shown are for the subsample corresponding 
to events with exactly one SVX-tagged jet. 
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Figure 7.2: Reconstructed-mass distributions from simulated it events for sev- 
eral input values for the top quark mass used in the simulation. The overlayed 
curves are predictions from the parametrization of templates at 18 different 
top mass values. The distributions shown are for the subsample corresponding 
to events with exactly two SVX-tagged jets. 




Figure 7.3: Reconstructed-mass distributions from simulated it events for sev- 
eral input values for the top quark mass used in the simulation. The overlayed 
curves are predictions from the parametrization of templates at 18 different 
top mass values. The distributions shown are for the subsample corresponding 
to events with one or more SLT-tagged jets and no SVX-tagged jets. 
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Figure 7.4: Reconstructed-mass distributions from simulated it events for sev- 
eral input values for the top quark mass used in the simulation. The overlayed 
curves are predictions from the parametrization of templates at 18 different 
top mass values. The distributions shown are for the subsample corresponding 
to events with no SVX-tagged or SLT-tagged jets, and the fourth jet having 
E T > 15 GeV and \r}\ < 2. 




Figure 7.5: Reconstructed-mass distribution for VT+jets events generated with 
the VECBOS Monte Carlo simulation. The smooth curves are the parametriza- 
tions of the reconstructed-mass distributions. The distributions are for SVX 
tagged, SLT tagged, and No Tag events used in the mass analysis. 
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Figure 7.6: The value of the top quark mass returned by the fitter for input 
top quark masses of 150, 175, and 200 GeV/c 2 . Each simulated experiment 
contains an admixture of signal and background events as described in the 
text. 
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Pull from Simulated Experiments 



Figure 7.7: The distribution of pulls between the returned value of the top 
quark mass and the input value of 175 GeV/c 2 . Each simulated experiment 
contains an admixture of signal and background events as described in the 
text. 
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Dashed line: Slope=1.0 (not a fit) 
Dotted line: Slope=0.92 ± 0.09 (fit) 
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Figure 7.8: The Gaussian width of the mass distribution as a function of the 
statistical uncertainty returned by the likelihood. Each simulated experiment 
contains an admixture of signal and background events as described in the 
text. The dashed line, which has a slope of 1.0, is not a fit. The dotted line is 
a fit to the points, and has a slope of 0.92±0.09. 
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RESULTS 



Having tested the mass likelihood procedure on simulated experi- 
ments, we now apply it to the data sample. Two separate fits are performed. 
The first is the standard mass likelihood fit defined in Section |7.2j . The second 
is the mass likelihood fit with the background fraction constraint removed. 
In this case, the background fraction is determined only from the shapes of 
the reconstructed-mass distributions for signal and background. After pre- 
senting these results, we check the consistency of the data with Monte Carlo 
expectation for ti plus background events in the expected proportion. 



8.1 Fits to data 

The mass likelihood procedure is applied to the events in the four 
mass subsamples. The measured values for the top quark mass for each sub- 
sample and the combined results are presented in Table [TT]. As discussed in 
Section [7.2| , the statistical uncertainties contain contributions from both the 
statistics in the data, the uncertainties in the expected background, and the 
uncertainties in the template fit parameters. The latter two contribute less 
than 1% to the total statistical uncertainty. Table |0 also shows the fitted 
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background fractions, which are constrained to the expected values via the 



background likelihoods in Fig. O. The mass fits for the four mass subsamples 
are statistically consistent with one another and are shown in Figure |S.1| . For 
each subsample, the background shape has been normalized to the fitted num- 
ber of background events via N% = x{ lt x N l ohs , and the signal plus background 
has been normalized to the number of data events {N l ohs ). The combined fit 



to all four subsamples is shown in Figure p72 . 

To investigate the impact of the background constraining term on 
the fitted top quark mass, we also performed mass likelihood fits with the 
constraint on the fraction of background removed. In this case, the shape of 
the mass distribution determines the background fraction. The results of the 



mass fits are presented in Table |8.2| . Several observations can be made from 
a comparison of these unconstrained mass fits with the constrained ones in 



Table |8J]. First, the tagged subsamples fit to zero background, although with 
large uncertainties, while the No Tag subsample yields a similar background 
content whether the background is free to float or not. Secondly, the masses 
show little sensitivity to removal of the background constraint. In general, 
one would expect the removal of the constraint to result in an increased sta- 
tistical uncertainty since information is being removed from the likelihood fit. 
For all subsamples however, the uncertainty in the mass decreases when the 
background constraint is removed. This is because the fitted number of signal 
events becomes larger. 

Since the background rates in the four mass subsamples are corre- 
lated, it is not correct to allow their background fractions to float relative to 
one another. On the other hand, it is reasonable to investigate whether the 
background constraint is affecting the top quark mass measurement. The re- 



sults in Tables |8.1| and |8.2| indicate very little sensitivity to the background 
constraint. 
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Data sample 


Number of events 


fit 

x b 


Top mass (GeV/c 2 ) 


SVX Double 


5 


0.03 ±0.02 


170.0l|j 


SVX Single 


15 


1 5 +a05 

u - id -0.04 


178.0+?;! 


SLT 


14 


0.53 ±0.09 


140.6±f2;| 


No Tags 


42 


0.693:3 


182.1^ 9 7 


All subsamples 


76 


0.51 ±0.06 


176.1+f; 2 , 



Table 8.1: Results of applying the mass likelihood procedure to the four sub- 
samples and for all subsamples combined. The background fractions are con- 
strained to their expected values via the curves in Fig. |4.1| . For the four 
subsamples, we show the fitted background fraction and the fitted top mass. 
Also shown is the final mass value obtained when combining all four subsam- 
ples. The combined background fraction is the average of the Xb fit results 
weighted by the number of events in the subsamples. 

8.2 Comparison of data to expectations 

Up to this point we have assumed that our data sample is a mixture 
of standard model ti signal plus background. Using simulated experiments 
(with M top =175 GeV/c 2 ), we quantify the probability that our data sample is 
consistent with a mixture of ti plus background with the background fractions 



given in Table |0 



We first check that the statistical uncertainty obtained from the data 
sample (5.1 GeV/c 2 ) is reasonable. Figure |Q| shows the distribution of statis- 



tical uncertainties from simulated experiments along with the value we obtain 
for our data sample. We find that 5% of simulated experiments yield a sta- 
tistical uncertainty equal to or smaller than the value from our data sample. 
While this number is small, it is reasonable. 

A further check is provided by comparing the minimum of the nega- 
tive log-likelihood obtained in the data sample with the values obtained from 
a large sample of simulated experiments (Fig. |S.4| ). A value of the negative 
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Data sample 


Number of events 


fit 

x b 


Top mass (GeV/c 2 ) 


SVX Double 


5 


n n +0 - 6 

u - u -0.0 


169.9l|? 


SVX Single 


15 




177.6t£? 


SLT 


14 


u - u -o.o 


146.2^ 


No Tags 


42 


0.53 ±0.22 


ISO.S^ 1 


All subsamples 


76 


0.29 ±0.20 


176.2 ±4.8 



Table 8.2: Results of applying the mass likelihood procedure to the four sub- 
samples and for all subsamples combined. The background fractions are free 
parameters in the mass likelihood fit. For the four subsamples, we show the 
fitted background fraction and the fitted top mass. Also shown is the fi- 
nal mass value obtained when combining all four subsamples. The combined 
background fraction is the average of the Xb fit results weighted by the number 
of events in the subsamples. 

log-likelihood larger than expected from simulated experiments might indicate 
that either the reconstructed-mass distribution is not well modeled or that the 
background fractions in the sample are not properly estimated. The distri- 
bution shows that the value obtained from the data is quite consistent with 
standard model tt plus background, as evidenced by the 79% probability of 
obtaining a value of — log £ larger than the one seen in the data. 



8.3 Results from 6-tagged events 

In previous publications [[[], ||, the top quark mass was measured 
using only events containing SVX and/or SLT tagged jets among the leading 
four jets as a single sample (7 events in Ref |J, 19 events in Ref @). The final 
sample of 34 6-tagged events has been analyzed as part of our four subsample 
fit using the likelihood method described in Sections |6| and [7]. The 34 tagged 
events may be treated as three subsamples or they may be fit as a single 34 
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event sample [S^j. The results of fitting the 34 6-tagged events as a single 
sample are shown in Fig. O. The likelihood fit yields a top quark mass of 
173.7±6.2(stat.) GeV/c 2 with a top fraction of 75%. Treating the 34 6-tagged 
events as three subsamples, we obtain a top quark mass of 174.0±5.7(stat.) 
GeV/c 2 . The 8% improvement is consistent with expectations from simulated 
experiments (see Section |3.3.1| ). 



8.4 Alternate techniques 

Other analyses which have different selection criteria and/or modi- 
fied formulations of the x 2 have been performed. The analyses are aimed at 
improving the probability for choosing the correct combination. The first of 
these analyses, the £** analysis, uses two additional terms in the x 2 to aid in 
choosing the correct combination. The second analysis uses a looser definition 
for 6-jet tagging to increase the number of double 6-tagged events. Values of 
the top quark mass from these two analyses are consistent with the results 
presented in this report and are summarized in Appendix [B|. 
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Figure 8.1: Results of applying the likelihood procedure to the four subsam- 
ples. The figure show the data (histogram), fitted background (shaded hatched 
region), and fitted signal (shaded non-hatched region). The insets show the 
shapes of — log C versus top mass, from which we extract the fitted top quark 
mass and its statistical uncertainty. 
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Figure 8.2: Result of applying the likelihood procedure to the combined sub- 
samples. The figure shows the data (histogram), fitted background (shaded 
hatched region), and fitted signal (shaded non-hatched region). The inset 
shows the shape of — log C versus top mass, from which we extract the best 
estimate of the top quark mass and its statistical uncertainty. 
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Figure 8.3: Distribution of statistical uncertainties from simulated experiments 
of ti (M top =175 GeV/c 2 ) plus background. Each simulated experiment con- 
tains a mixture of signal and background events as described in the text. Also 
shown is the statistical uncertainty obtained from our data sample. The prob- 
ability for obtaining a smaller uncertainty in the simulated experiments is 5%. 
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Figure 8.4: Distribution of the minimum value of — log C from simulated exper- 
iments of ti (M top =175 GeV/c 2 ) plus background. Each simulated experiment 
contains a mixture of signal and background events as described in the text. 
Also shown is the minimum — log£ value obtained from our data sample. The 
probability for obtaining a larger value of — log C in the simulated experiments 
is 79%. 
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Figure 8.5: Result of applying the likelihood procedure to the 34 6-tagged 
events, treated as a single sample. The figure shows the data (histogram), 
fitted background (shaded hatched region), and fitted signal (shaded non- 
hatched region). The inset shows the shape of — log£ versus top mass, from 
which we extract the top quark mass and its statistical uncertainty. 



Section 9 

SYSTEMATIC 
UNCERTAINTIES 

The general procedure for estimating the systematic uncertainty on 
M top from a given source is handled similarly for all sources. For a given source 
of uncertainty, we make a change in the input value, and evaluate the impact 
on the measured top quark mass. The change is either a one standard devi- 
ation (la) uncertainty on the variable in question, or a change in an input 
assumption. The change in the top quark mass is evaluated using simulated 
experiments (see Section |7.^ ). We perform a large number of simulated ex- 
periments with (a) the nominal input value of the variable or the standard 
assumption, and (b) a "la" shift in the variable value or the changed assump- 
tion. The reconstructed-mass distribution from each simulated experiment is 
fit using the same likelihood procedure as used on the data sample, thus ob- 
taining a measured top quark mass. The likelihood procedure includes the 
same templates as used with the data. The systematic uncertainty is defined 
as the difference in the median top quark mass between the two ensembles (a) 
and (b). The distribution of reconstructed top quark masses from simulated 
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experiments in which all inputs are set to their nominal values was shown in 
Fig. 0. 

9.1 Uncertainties from the energy scale 

The event reconstruction algorithm varies the measured momenta of 
the jets, leptons, and unclustered energy to fit the kinematics of the hypoth- 
esized ti decay. The energy scale for electrons and the momentum scale for 
muons are known to a precision of 0.14% and 0.065% respectively [plj] . This 
uncertainty has a negligible effect on the uncertainty in the top quark mass 
measurement. The uncertainty on the quantity Xt, the transverse energy 
beyond the partons associated with the ti event, has been discussed in Sec- 
tion f>.2[ This uncertainty is large, but large variations of Xt do not have a 
significant impact on the mass reconstruction. The % T is evaluated through 
the measurements of the charged lepton, the jets, and the unclustered energy, 
and is therefore not an independent measurement (see Section 5.2). To avoid 



correlations it is not used as a measurement in the fit, but is used to estimate 
a starting value for the transverse momentum of the neutrino. Thus, the en- 
ergy scale uncertainty in the measured top quark mass is dominated by the 
uncertainty in the measurement of the jet momenta. 

The total uncertainty in the jet Pt scale is taken as the quadrature 



sum of all uncertainties discussed in Section |5.1| . We apply +la and —la shifts 
to the jet momenta in ti signal and background events, and measure the effect 
on the measurement of the top quark mass. For the SVX Single subsample, 
the distributions of reconstructed masses for —la and +la shifts in the Pt 
scale are shown in Figs. (9.11 (a) and (c) respectively. These distributions may 
be compared to Fig. |9.1j (b) which shows the distribution obtained from the 
default momentum scale. As expected, a clear shift in the reconstructed-mass 
spectrum is observed. We generate analogous distributions for the other three 
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mass subsamples and for the background mass distribution. To obtain the sys- 
tematic uncertainty, we generate two large samples of simulated experiments. 
In the first sample, we choose the reconstructed masses for tt events at random 
from distributions like the one in Fig. |9j](a). In the second sample, we use 
distributions like the one in Fig. |9.1| (c). The simulated experiments in each of 
these samples are fit using the standard templates and the likelihood technique 
described in Section [7|. The median top mass from the simulated experiments 
in the two samples differ because of the applied jet Pt scale shifts. The distri- 
bution of reconstructed top masses from the two (jet Pt shifted) samples are 



displayed in Fig. |9]2|. We take half the difference between the medians of the 
— la and +la distributions (from Fig. |9.2| ) as the uncertainty on the top quark 
mass measurement due to the Pt scale uncertainty. Using this prescription, 
we obtain a top mass uncertainty of ±4.4 GeV/c 2 from the jet Pt scale. 



9.2 Initial and final state hard radiation 

QCD radiation that produces jets can originate from the outgoing 
(final state) partons, the incoming (initial state) partons, or from interference 
among the two. The interference effect is expected to be small [SB] and is not 
considered here. 

The effects of initial state radiation (ISR) and final state radiation 
(FSR) on the measurement of M top are studied using the PYTHIA program 
since it allows the two effects to be studied in isolation from one another. The 
approach used to evaluate the systematic uncertainty due to ISR is to compare 
the median mass from simulated experiments using the standard PYTHIA 
settings to the median mass from simulated experiments with ISR turned off. 
The median mass from simulated experiments for the no-ISR PYTHIA sample 
is found to be lower than that of the standard PYTHIA sample by 2.6 GeV/c 2 . 
The uncertainty is taken to be one half of the shift in median mass between 
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Figure 9.1: Distribution of reconstructed masses for SVX Single tagged events 
from the HERWIG tt Monte Carlo simulation using an input top mass 175 
GeV/c 2 for (a) a —la shift in the jet Pt scale, (b) no shift in the jet Pt scale, 
and (c) a +lcr shift in the jet Pt scale. These distributions are used as inputs 
to generate the samples of simulated experiments described in the text. 



106 



1000 



800 



Jet P T Shift 



Median (GeV/c) = 171.2 



+ 

180 



> 

<u600 

o 

CN 

B 

«400 



200 



150 



160 



170 



180 



190 



200 



Fitted Top Mass from Simulated Experiments (GeV/c ) 



Figure 9.2: Distribution of fitted top masses returned from the likelihood pro- 
cedure for simulated experiments. The solid histogram shows the distribution 
when the jet Pt scale is shifted down by —la, and the dashed histogram 
shows the results when the jet Pt scale is shifted up by +la. The median top 
masses for each are indicated on the figure, from which we obtain a systematic 
uncertainty of 4.4 GeV/c 2 . 
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the standard PYTHIA simulation and the no-ISR PYTHIA simulation, which 
is 1.3 GeV/c 2 . We assume that the shift is symmetric with the amount of ISR 
and therefore the uncertainty is ±1.3 GeV/c 2 . 

Extracting the effects due to final state radiation is a more subtle 
exercise because PYTHIA, like HERWIG, describes jet formation through a 
parton shower. The effects of modeling the softer components on the mea- 
surement of M top are described by the studies of soft gluon radiation (see 
Section |5.1.1| ). In this discussion, we are therefore referring to the "harder" 
component of FSR, which leads to extra jets in the final state. To isolate 
the effect of FSR, we use a sample of PYTHIA events which have ISR turned 
off. We select a subsample of these events that have exactly four jets ( either 
four high-E T jets, E T > 15 GeV and \r)\ < 2.0, or three high-E T jets plus 
one with E T > 8 GeV and \rj\ < 2.4), all of which can be uniquely matched 
to the partons from ti decay (within a cone of 0.4). Using the procedure of 
simulated experiments, we take the systematic uncertainty to be half the dif- 
ference between the no-ISR PYTHIA events with no restriction on the number 
of jets and the subsample of events with exactly four jets uniquely matched 
to the partons from ti decay. We assume that this difference is approximately 
independent of the amount of ISR present. The median top mass from the 
no-ISR sample with exactly four jets is found to be larger than the default 
no-ISR sample by 4.4 GeV/c 2 . The systematic uncertainty is taken to be one 
half of the difference in the median top masses between the two samples, or 2.2 
GeV/c 2 . As with ISR, we assume that the shift in top mass is symmetric with 
respect to the amount of FSR, so the systematic uncertainty is ±2.2 GeV/c 2 . 



9.3 Background mass distribution 



In generating the default background distributions with the VEC- 
BOS program, we used the VF±3 partons matrix elements and chose a scale of 
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Q 2 = (Pt) 2 - This Q 2 scale is not only used in VECBOS for the computations 
of the matrix elements and the evaluation of the parton distribution functions, 
but it is also used in the evolution of the parton showers to limit the Pt of 



additional jets |p9| . As a result, the shape of the reconstructed-mass distribu- 
tion is sensitive to the choice of scale. The systematic uncertainty from using 
the assumed background shape was evaluated by changing the Q 2 scale from 
(Pt) 2 to Myy- Simulated experiments using ti and the modified background 
shape (Q 2 = Myy) were fit to the default signal and background probability 
distributions as described in Section [7]. The median mass from these simu- 
lated experiments was found to differ by 1.3 GeV/c 2 from simulated experi- 
ments using the default background shape. The systematic uncertainty from 
the background shape modeling is taken to be this difference symmetrized, or 
±1.3 GeV/c 2 . 



9.4 b- Tagging 

A systematic uncertainty in the top mass measurement may arise 
from an uncertainty in the SVX and SLT tagging efficiencies. For SVX tagging, 
the primary uncertainty comes from the possible Et dependence of the SVX 
tagging efficiency which may differ from the simulation. The SVX tagging 
efficiency in data relative to the simulation is parametrized as a function of 
the jet Et, and is nearly flat. The systematic uncertainty is evaluated by 
assuming the largest possible Et dependence of this ratio given the available 
statistics. Comparison of results obtained using a flat ratio to those obtained 
with a sloped one gives a mass shift of only 0.1 GeV/c 2 

The jet Et dependence of the SLT tagging efficiency is better known 
than in the SVX case. However, a systematic uncertainty does arise from the 
uncertainty in the expected ratio of true to fake SLT tags in ti events. To 
estimate the sensitivity of our top quark mass measurement to this ratio, we 
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generate Monte Carlo it events in which all SLT tags are either (a) true tags, 
or (b) fake tags. We then produce two large samples of simulated experi- 
ments, each having the same number of observed events as in our four mass 
subsamples and including the appropriate background contributions. The two 
samples of simulated experiments consist of either events all from set (a) or 
all from set (b). The median top quark mass values from the two samples of 
simulated experiments differ by 0.8 GeV/c 2 . We take half of this difference, 
or 0.4 GeV/c 2 , as the corresponding systematic uncertainty in the top quark 
mass. Combining the systematic uncertainties from SVX and SLT tagging, we 
find an overall systematic uncertainty of ±0.4 GeV/c 2 . 



9.5 Parton distribution functions 

All of the Monte Carlo samples used to measure the top mass were 
generated with the MRSDO' [37|| set of parton distribution functions (PDF). 



This was the preferred PDF at the time the samples were generated. Newer 
distribution functions now exist, in particular ones which fit CDF's inclusive jet 
cross section. One such PDF, CTEQ4L provides a higher gluon content 
at lower momentum fraction than MRSDO'. We take the shift in the median 
top mass between samples generated with the two PDF's as the relevant un- 
certainty. We therefore assign a systematic uncertainty of ±0.3 GeV/c 2 in the 
top quark mass from this source. 



9.6 Monte Carlo generators 

The effect of using different Monte Carlo generators has also been 
studied. Previously, this was evaluated from the difference between the HER- 
WIG and ISAJET simulations. Because of the evidence that independent frag- 
mentation does not reproduce some aspects of the data (energy flow around 
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and between jets [ |6"T| , |!|) we will not use IS A JET here. We evaluate the sys- 
tematic uncertainty from the choice of Monte Carlo generators via the mass 
shift between the HERWIG and PYTHIA simulations. This gives a systematic 
uncertainty of ±0.1 GeV/c 2 . 



9.7 Summary of systematic uncertainties 

The relevant systematic uncertainties studied for the top mass mea- 
surement are listed in Table |9.1| . Combining all of these effects in quadrature 
gives a total systematic uncertainty of ±5.3 GeV/c 2 , or ±3% of 176.1 GeV/c 2 . 



Source 


Uncertainty (GeV/c 2 ) 


Jet energy measurement 


4.4 


Initial and final state radiation 


2.6 


Shape of background spectrum 


1.3 


b- Tagging 


0.4 


Parton distribution functions 


0.3 


Monte Carlo generators 


0.1 


Total 


5.3 



Table 9.1: Systematic uncertainties on the measurement of the top quark mass 
for this analysis. 



Section 10 



COMBINED TOP QUARK 

MASS 

The most precise measurement of the top quark mass in any single 
decay channel is obtained with events in the lepton+jets topology. The anal- 
ysis of such events leads to a mass of 176.1 ± 5.1(stat.) ± 5.3(syst.) GeV/c 2 . 
Measurements in the all-hadronic |T(| and dilepton |62[ decay topologies have 
also been made and can be combined with the lepton+jets result to reduce 
the overall uncertainty. Here we make some brief remarks on these analyses, 
and describe how the three measurements were combined. 

10.1 All-hadronic topology 

The top quark mass measurement in the all-hadronic topology used 
a sample of 136 events that satisfied several selection criteria, including the 
requirement of six or more jets, at least one of which was tagged as a b by the 
SVX. The estimated background in the sample was 108±9 events. The method 
for extracting a top mass was similar to the one used for the lepton+jets 
topology, and included a kinematic fit to each event and a likelihood fit to 



111 



112 



the resulting reconstructed-mass distribution. The results of the likelihood fit 
yielded a measured top quark mass of 186 ± 10(stat.) ± 12(syst.) GeV/c 2 [ITU 



A reevaluation of the systematic uncertainty on this measurement has led to a 
more accurate estimate of 5.7 GeV/c 2 f62fl . Appendix describes the details 
of this reevaluation. 

10.2 Dilepton topology 

The dilepton topology includes ti events in which the W + and W~ 
bosons each decay into an ev or [iv final state. The presence of two neutrinos, 
which are not observed in our detector, prevents a straightforward event-by- 
event kinematic fit to the ti decay hypothesis. Therefore, we have measured the 
top quark mass from dilepton events using a weighting method [13|, |64|. In 
this method the vector sum of the neutrino transverse momenta, as predicted 
after making certain assumptions, is compared to the observed missing trans- 
verse momentum [|6^| . From a sample of eight events with an estimated back- 
ground of 1.3±0. 3 events we obtain a mass of 167.4 ± 10.3(stat.) ± 4.8(syst.) GeV/c 2 
A brief description of the method, and some additional information not re- 



ported in Ref. p2J is given in Appendix O. 

10.3 Combining the measurements 

Each of the three top quark mass measurements is associated with 
a statistical and systematic uncertainty. The statistical uncertainties are un- 
correlated, since the samples are statistically independent. However, the sys- 
tematic uncertainties are correlated, and these correlations must be included 
when combining the results. 

The systematic uncertainties in the measurements from each decay 
topology |62|] are assigned to one of five independent categories: 



113 



1. jet energy scale, 

2. signal model (ISR, FSR, PDF, 6-tagging), 

3. Monte Carlo generator, 

4. background model, 

5. Monte Carlo statistics. 

The assignment of the systematic uncertainties for each of the three mass 



analyses to these categories is shown in Table |1(J.1| . In the lepton+jets mea- 
surement, the statistical uncertainty in the Monte Carlo simulation is included 
in the global statistical uncertainty 





Systematic uncertainty (GeV/c 2 ) 


Systematic category 


Lepton+jets 


All-hadronic 


Dilepton 


Jet energy scale 


4.4 


5.0 


3.8 


Signal model 


2.6 


1.8 


2.8 


M.C. generators 


0.1 


0.8 


0.6 


Background model 


1.3 


1.7 


0.3 


M.C. statistics 




0.6 


0.7 


Total 


5.3 


5.7 


4.8 



Table 10.1: Systematic uncertainties for each of the three mass analyses 
grouped into the five categories. Also shown is the total systematic uncer- 
tainty for each analysis. 

For each of the five categories, the systematic uncertainties in each 
of the three measurements are assumed to be either uncorrelated or 100% 
correlated. The jet energy scale uncertainty is taken to be 100% correlated 
since all three analyses use the same detector and the same jet clustering al- 
gorithm. The systematic uncertainties coming from the signal model and the 
Monte Carlo generator are also assumed to be 100% correlated since all three 
analyses use the HERWIG Monte Carlo generator to simulate ti events. The 
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uncertainties in the background shape are assumed to be uncorrected because 
the background processes for each analysis are different. The correlation coef- 
ficients between the three pairs of analyses are given in Table |10.2| . 



Correlation coefficients 



Systematic category LJ/AH LJ/LL AH/LL 



Jet energy scale 


1.0 


1.0 


1.0 


Signal model 


1.0 


1.0 


1.0 


M.C. generators 


1.0 


1.0 


1.0 


Background model 


0.0 


0.0 


0.0 


M.C. statistics 


0.0 


0.0 


0.0 



Table 10.2: Correlation coefficients between the three mass analyses for the 
five categories of systematic uncertainty. Here, LJ signifies the lepton+jets 
analysis, AH the all-hadronic analysis, and LL the dilepton analysis. 

The inputs into the calculation for combining the mass measurements 
are the three top quark mass measurements cited in this section, their statis- 
tical uncertainties, and the systematic uncertainties and their correlations as 



listed in Tables |10.1| and |10.2| . The calculation uses a generalized chisquare 



method with full covariance matrix (see for example ||65||), and yields: 

m t = 176.1 ±6.6 GeV/c 2 . (10.1) 

In the calculation, the central value can be written as the weighted sum of 
the three input central values. The weights, which depend on the statistical 
and systematic uncertainties and the correlations, are found to be 0.65 (lep- 
ton+jets), 0.19 (dilepton), and 0.16 (all-hadronic). If we define a statistical 
uncertainty on the combined result as the quadrature sum of the weighted 
individual statistical uncertainties, that combined statistical uncertainty is 
±4.2 GeV/c 2 . The combined systematic uncertainty, defined as the quadrature 
difference between the total and statistical uncertainties, is then ±5.1 GeV/c 2 . 
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SUMMARY 



The first evidence of the production of top quark pairs in pp colli- 
sions was reported by CDF in 1994 [|IJ. From a sample of seven candidate 
lepton+jets events with an expected background of 1.4 events, the top quark 
mass was measured to be 174±10(stat)ti2(syst) GeV/c 2 . Since that time, 
both CDF and DO have analyzed their full Run 1 data samples and have pub- 
lished their measurements of the ti production cross section |7|, and top 
quark mass []]. Indirect measurements of the top quark mass using data 
from LEP and SLC have been made fllT] , and are consistent with the direct 
measurements, although with a substantially larger uncertainty. 

This report has described in detail the best single channel measure- 
ment of the top quark mass. A letter on the measurement is already pub- 
lished ||. The likelihood method uses parametrized templates, which results 
in a continuous likelihood shape as a function of the top quark mass from which 
the top quark mass and statistical uncertainty are evaluated. The statistical 
precision of the top quark mass measurement has benefitted from a larger 
data sample than earlier measurements Jl|, [|, and through subdivision of the 
data sample into non-overlapping subsamples according to the 6-tagging infor- 
mation. Systematic uncertainties have been considerably reduced, primarily 
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through a better understanding of the jet energy measurements which resulted 
in smaller uncertainties on jet Pt- 

From the 106 pb _1 Run 1 data sample, we measure the top quark 
mass in the lepton+jets topology to be 176.1±5.1(stat.)±5.3(syst.) GeV/c 2 . 
Measurements of the top quark mass in the all-hadronic [[HJ and dilepton p2[ 
decay topologies are consistent with this measurement. Combination of the 
three measurements from CDF gives a top quark mass of 176.1±6.6 GeV/c 2 . 
The DO collaboration has also published results on the top quark mass mea- 
surement in the lepton+jets and dilepton channels, from which they obtain a 
combined top quark mass of 172.1 ±7.1 GeV/c 2 p6| . The measurements of the 
top quark mass from the CDF and DO experiments are consistent with each 
other, therefore, their Run 1 measurements have been combined to obtain a 
top quark mass at the Tevatron of 174.3 ± 5.1 GeV/c 2 |67]]. This measurement 
represents the most precise measurement of any of the quark masses. 
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Appendix A 

EVENTS IN THE MASS 
ANALYSIS 



The individual reconstructed masses of all events in the four subsam- 
ples are listed in Tables |A.1| through |A.4. 
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Mass 

Run Event (GeV/c 2 ) 

40758 44414 175.3 

67824 281883 170.1 

65581 322592 152.7 

67971 55023 183.5 

68464 547303 151.1 



Table A.l: List of events which are in the SVX Double subsample. Shown are 
the run and event numbers and the reconstructed top mass for the solution 
having the lowest \ 2 ■ 



Mass 



Run 


Event 


(GeV/c 2 ) 


43096 


47223 


288.6 


45610 


139604 


180.0 


45879* 


123158 


180.1 


59698* 


31639 


187.4 


63247 


65096 


161.0 


63641 


3054 


173.3 


68006 


44672 


243.4 


64901 


569801 


156.3 


69683* 


135095 


163.2 


56911* 


114159 


156.7 


67515 


298909 


174.6 


68312* 


821014 


202.4 


68739 


425355 


170.9 


69781* 


266905 


182.8 


56669 


21631 


152.1 



Table A. 2: List of events which are in the SVX Single subsample. Shown are 
the run and event numbers and the reconstructed top mass for the solution 
having the lowest X Events labelled with a * have both SVX and SLT tagged 
jets. 
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Mass 



Run 


Event 


(GeV/c 2 ) 


45705 


54765 


186.3 


45880 


31838 


130.4 


43351 


266423 


162.4 


66368 


91765 


137.9 


66500 


421896 


173.0 


67879 


30394 


141.1 


69005 


181134 


129.6 


58908 


41102 


138.6 


60998 


423792 


162.0 


61334 


57897 


183.1 


64721* 


229200 


181.0 


65298* 


747402 


149.4 


65648 


203840 


203.2 


67515 


616477 


149.9 



Table A. 3: List of events which are in the SLT subsample. Shown are the run 
and event numbers and the reconstructed top mass for the solution having the 
lowest x 2 - Events labelled with a * have two SLT-tagged jets. 
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Mass 






Mass 


Run 


Event 


(GeV/c 2 ) 


Run 


Event 


(GeV/c 2 ) 


46492 


57501 


179.2 


58696 


83095 


137.6 


41301 


45902 


175.7 


59948 


105232 


115.4 


43421 


65648 


147.8 


60634 


350037 


151.2 


47757 


262594 


219.6 


61167 


332223 


167.3 


45757 


30003 


173.0 


63265 


5385 


255.2 


45144 


107403 


189.2 


64041 


473567 


247.5 


60656 


96710 


180.3 


64997 


78806 


192.0 


60746 


121257 


180.1 


65179 


215794 


195.7 


61511 


75858 


113.0 


67391 


50780 


184.9 


62981 


85084 


125.0 


67757 


631972 


172.0 


64861 


121618 


178.8 


68144 


100373 


178.3 


64934 


400688 


215.4 


68231 


78554 


177.7 


66046 


507038 


164.2 


68374 


312573 


139.1 


66207 


12039 


154.4 


68553 


707057 


130.4 


66315 


365275 


230.3 


68570 


897728 


142.6 


67862 


631243 


114.2 


68593 


88427 


144.0 


68006 


176291 


120.9 


69519 


430034 


160.0 


68939 


352425 


173.1 


70000 


26023 


161.1 


69520 


307639 


235.2 


57438 


71994 


253.1 


70578 


351956 


143.0 


64901 


505659 


108.1 


70986 


227609 


176.2 


67397 


105755 


190.0 



Table A. 4: List of events which are in the No Tag subsample. Shown are 
the run and event numbers, and the reconstructed top mass for the solution 
having the lowest x 2 . 



Appendix B 



ALTERNATE MASS 
ANALYSES 

A number of alternate mass analyses have been performed using the 
Run 1 data sample. We discuss two alternate analyses which are aimed at 
improving the statistical and/or systematic uncertainty on the top quark mass 
measurement using some subsample of events. Another goal is to check our 
default technique by employing complementary strategies by: (i) using more 
event information associated with b tagging and jet charge, and (ii) reducing 
the sample to the most complete events, ie. those where we have two b-tagged 
jets. The first of these techniques includes additional terms in the likelihood 
function, which improves the probability for choosing the correct jet-to-parton 
configuration at the expense of reduced statistics. The second technique uses 
three 6-tagging algorithms to explore a subsample of the data set that consists 
of events with two 6-tagged jets among the leading four jets. Neither of these 
two techniques is found to yield a more precise measurement than the mass 
analysis described in the body of this report. In this appendix, we briefly 
describe these two mass analyses. 
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B.l The C** fitting technique 



The £** technique |k| aims at improving the fraction of correct jet- 
to-parton assignments by combining three independent sources of event infor- 
mation into a single parameter. These sources are: 

• x 2 f° r it-like kinematics as described in Section || 

• probability for the jets assigned as b jets to originate from b quarks, and 
the two jets assigned to the hadronic W-decay to originate from light 
quarks. The probability is evaluated using the jet probability (JPB) 
algorithm 



probability to observe a given jet charge [[nj for b and b quarks in tt 
events. 



B.l.l Definition of 

The JPB algorithm evaluates for each charged track in a jet the 
probability that it comes from the primary vertex. For each jet the track 
probabilities are combined into an overall probability (JPB) that the jet is 
consistent with the zero lifetime hypothesis. Due to the long lifetime of b 
hadrons, the JPB distribution for b quark jets exhibits a strong peak near 
zero. Non-6 jets in tt events are produced either through the decays of W- 
bosons to (u, d) and (c, s) quark pairs, or production of gluon jets from initial 
or final state radiation. With the exception of the charmed quarks, the non-6 
jets exhibit a flat JPB distribution. The charm quark jets produce a small 
peak near zero which can be ignored given its relative size. Unless otherwise 
noted, charm quark jets are understood to be included in the "non-6" quark 
distribution of JPB. 

We incorporate the JPB variable into the x 2 definition by introducing 
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the following selection function: 

C = x 2 - 2 • In [P(JPBi) • P(JPB 2 ) • P(JPB 3 ) • P(JPB 4 )] . (B.l) 

The x 2 is the same as the one defined in Section |6], and V(JPB[) is the prob- 
ability density for the i th -jet assignment (i=l,...4). The V functions in £* 
depend only on jet type, since one function is appropriate for both b and b 
jets, and another for non-6 quarks. While the x 2 value is in general different 
for each of the 24 combinations, only six distinct values occur for the second 
term in C*. Groups of four, corresponding to the interchange of the b and b 
quarks (and the two neutrino Pz solutions), have the same contribution from 
this second term. 

We used the HERWIG Monte Carlo and the full CDF detector sim- 
ulation to generate the (b and b) and non-(6 and b) f72| probability density 
distributions. We only considered events in which the leading four jets corre- 
sponded to the four primary partons from ti decay, which limits us to 56% of 
the sample. Of this subset, we found that the largest fraction of correct as- 
signments based on selecting combinations with minimum £* was 48%, which 
was obtained with a jet clustering cone size of 0.4 and a minimal track Pt of 
1.0 GeV/c. 

To incorporate additional information pertaining to the charge of the 
b and b jets, we define a new selection function, 

£** = £*- 2- In [C(Qb)-C(Q 6 )], (B.2) 

where C(Qj et ) is the jet charge probability density. The jet charge is defined 
as in Ref. 0: 

Q jet = E£? g< - IP. ■ el* (R3) 

where e is the unit vector along the jet axis, qt and pi are the charge and 
momentum of the i th track, and the sum extends over all n tr ^ charged particles 
in a fixed cone around the jet. To determine optimal choices for the cone size 
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and the weighting factor k, we varied the jet cone size from 0.35 to 1.0 and 
k from 0.4 to 1.2, and compared the significance of separation between the 
b and b C(Qj et ) distributions. The results were relatively insensitive to the 
exact values of these parameters. Since we found no strong dependence on 
these parameters, we chose the same cone size as used to calculate the JPB 
probability and for simplicity selected k — 1. 

B.1.2 Event selection and number of expected back- 
ground events 

In this analysis we select events with at least one SVX or SLT tag. 
All of the standard lepton and jet corrections discussed previously in this paper 
are applied. A total of 34 events are accepted which are identical to the tagged 
events shown in Table |3.1| . Since the JPB algorithm uses tracks reconstructed 
in the silicon vertex detector, we require that each event has at least one jet 
with associated SVX tracks. We also require that the combination with the 
lowest value of £** has a x 2 ( as defined in Section |D less than 10. Only 
solutions in which a 6-tagged jet is assigned to a b parton are considered. We 
find that 27 of the 34 events pass these requirements. 

We take the combination with the lowest £** value as the most likely 
decay chain of the ti into the four highest Et jets. Monte Carlo studies show 
that switching from x 2 to the £** selection increases the probability of making 
the correct jet-to-parton assignments. The probability of correctly assigning 
the four highest E T jets to the ti daughter partons increases from 30.5±0.7% 
to 37.3±0.6%. This fraction is "a priori" limited to a maximum of 56%, due 
to jets from ISR and FSR. 

The number of expected background events for the 34 tagged events 
is estimated to be 10.2 ± 1.5, which includes a background of 7.6 ± 1.3 for 
the 14 events with only SLT tags. This analysis reduces the number of SLT 
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tagged events from 14 to 7 (no SVX tagged events are cut out). Using the 
method described in Section |47|, the expected background for the 7-event SLT 
sample was evaluated to be 3.2tol events. We therefore calculate an expected 
background for the 27-event sample of 5.8±o;g events, which corresponds to a 
background fraction Xb =0.21^o;q3. 



B.1.3 Result of the likelihood fit 

The evaluation of the top mass uses the same techniques described in 
Section ^. The result of the fit is shown in Fig. |B.1| . The histogram represents 



the reconstructed mass distribution for the 27 data events. The shaded area 
corresponds to the background fraction returned by the fitting procedure, and 
the smooth curve shows the sum of the fitted background and signal contribu- 
tions. The insert displays the likelihood shape with the background fraction 
constrained to 0.2lio;o3. The resulting fit yields: 

M top = 170.3t^(stat.) GeV/c 2 . (B.4) 

The soundness of the procedure was tested using simulated experi- 
ments. Figure |B.2| (a) shows the pull distributions for simulated experiments, 
and Fig. [B.2| (b) shows the average (of the positive and negative) statistical 



uncertainty returned from the likelihood fit. The arrow indicates the fit result 
from the data sample. We find that 44% of simulated experiments have a 
statistical uncertainty smaller than measured in the data sample. 

Using simulated experiments, we compared the expected statistical 
uncertainty from 34 tagged events using the standard kinematic fit with 27- 
event experiments using the £** technique. The studies indicated that for 
samples of this size, we could reduce the top quark mass measurement uncer- 
tainty by ~ 0.5 GeV/c 2 over the standard kinematic \ 2 i if we consider the 34 
events as a single sample. 
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B.1.4 Systematic uncertainties 

The same categories of systematic uncertainties which were shown in 
Section ^| are present in this analysis. Moreover, we introduce a new systematic 
uncertainty which accounts for a possible difference in the tracking efficiency 
between data and simulation. This uncertainty is introduced because both 
the JPB algorithm and the jet charge calculation have some sensitivity to 
the tracking efficiency in jets. Using simulated experiments, we find an ex- 
pected uncertainty in the top quark mass of 0.9 GeV/c 2 from this source. The 



systematic uncertainties for the £** method are summarized in Table B.l 



Systematic 


Value 


uncertainty 


(GeV/c 2 ) 


Jet energy measurement 


4.0 


Initial and final state radiation 


2.7 


Shape of background spectrum 


0.5 


b- Tagging 


0.3 


Parton distribution functions 


0.6 


Monte Carlo generators 


0.8 


Tracking efficiency 


0.9 


Total 


5.1 



Table B.l: Systematic uncertainties for the C ** analysis. 



In conclusion, the £** analysis technique has been applied to a 27- 
event subset of the 34 tagged events, and leads to a top quark mass mea- 
surement of 170. 31^4 (stat.) ± 5.1 (syst.) GeV/c 2 . This value is in good 
agreement with the results presented in Section |8|. 
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100 150 200 250 300 350 

Reconstructed Mass (GeV/c 2 ) 

Figure B.l: Results of the £** likelihood fit to 27 6-tagged events (histogram). 
The shaded area corresponds to the background returned by the fit (5.2 ± 1.4) 
and the smooth line is the sum of fitted signal and background contributions. 
The inset shows the shape of the -log £** versus top mass from which we 
extract the fitted top quark mass and background contribution. 
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-3-2-10123 

Pull from Simulated Experiments 




GeV/c 



Statistical Uncertainty from Simulated Experiments 

Figure B.2: Results from 500 simulated experiments using the £** likelihood 
technique. Each simulated experiment consists of 27 events, as in the data 
sample. The upper figure shows the pull distribution between the fitted top 
mass and the input value (175 GeV/c 2 ), and the bottom figure shows the aver- 
age statistical uncertainty on the fitted top mass as returned by the likelihood 
procedure. The arrow represents the measured value in the data sample. 
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B.2 Fitting double 6-tagged events 



This analysis considers only events which contain two 6-tagged jets [173] 
To increase the acceptance for double 6-tagged events beyond what the SVX 
and SLT algorithms provide, we allow one of the b jets to be tagged by the 
JPB algorithm. Because the JPB variable depends upon the impact param- 
eters of the tracks in the jet with respect to the primary interaction vertex, 
this algorithm is correlated with the SVX tagging algorithm. We expect little 
or no correlation between the JPB and SLT tagging algorithms. We consider 
a jet 6-tagged if it has a JPB value less than 5%. 

Events are selected using the same selection criteria described in Sec- 



tion y. After we apply analysis cuts 1-7 from Section |3.3| and require that 
two jets are tagged by at least one of the three 6-tagging algorithms, the data 
sample consists of 11 events. Monte Carlo simulations showed that the res- 
olution on the measured top quark mass can be improved by requiring the 
invariant mass of the two untagged jets to be near the W^-boson mass. A cut 
of 60 < Mjj < 100 GeV/c 2 was found to yield the lowest uncertainty on the 
measured top quark mass. Nine of the eleven events are found to survive the 
W mass cut. The nine events are a subset of the sample of 34 tagged events. 
The reconstructed top masses of these nine events are listed in Table |B.2| and 
are shown in Figure B.3| . 

The expected backgrounds are estimated to be 0.22±0.08 events from 
Wbb + Wcc processes, 0.05 ±0.02 events from non-W background (e.g. bb 
production), and 0.13 ±0.05 events for non-heavy flavor background such as 
WW and WZ processes. The total number of background events is then 
estimated to be 0.4±0.1. 

The method for evaluating the top quark mass from this data sample 
is the same as the procedure discussed in Section |7|. The results of the fit are 
shown in Fig. [B.3| . The figure shows the mass distribution obtained from data 
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overlayed with the fitted results from the Monte Carlo simulation. The inset 
shows the distribution of —A log £ as a function of the top mass for the nine 
data events. 

The evaluation of systematic uncertainties are carried out in a similar 
manner to that which was discussed in Section ^|. The results are shown in 
Table |B.3| . The uncertainty due to background shape is appreciably reduced 
compared to the four subsample analysis because of the smaller background 
fraction. 

Using the techniques described in this section on the nine double 
tagged events, we measure the top quark mass to be 171.8±7.2(stat.)±4.3(syst.) 
GeV/c 2 . This measurement is consistent with the results presented in Sec- 
tion M. 



Di-jet mass Top mass 



Run 


Event 


Tags 


(GeV/c 2 ) 


(GeV/c 2 ) 


40758 


44414 


svx + svx 


83.9 


175.4 


59698 


31639 


SVX + (SLT & JPB) 


79.5 


187.4 


63247 


65096 


SVX + JPB 


81.3 


161.0 


64721 


229200 


SLT + SLT 


81.6 


181.0 


65298 


747402 


SLT + JPB 


60.0 


149.4 


65581 


322592 


(SVX & SLT) + SVX 


66.2 


152.7 


67824 


281883 


(SVX & SLT) + SVX 


73.3 


170.1 


67971 


55023 


SVX + SVX 


98.1 


183.5 


68464 


547303 


SVX + svx 


87.3 


151.1 



Table B.2: List of events used in the double 6-tagged analysis. Shown are 
the run-event numbers, the algorithms which tagged the two jets, the di-jet 
mass of the two untagged jets, and the reconstructed top mass for the solution 
having the lowest \ 2 . If a jet is tagged by two different algorithms, both tags 
appear in parentheses. 
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Systematic uncertainties 


Values (GeV/c 2 ) 


Jet energy measurement 


4.1 


Initial and final state radiation 


1.1 


Shape of background spectrum 


<0.1 


b- Tagging 


0.4 


Parton distribution functions 


0.3 


Total 


4.3 



Table B.3: Summary of systematic uncertainties in the top mass measurement 
from double 6-tagged events. 
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Figure B.3: Distribution of reconstructed mass for the nine data events and 
Monte Carlo simulation for the Double b-tag analysis. The background dis- 
tribution (hatched) has been normalized to the expected background of 0.4 
events. The sum of signal+background (dashed line) has been normalized to 
the fitted number of tt and background. The inset shows the shape of the 
likelihood function versus the top quark mass, from which we extract the top 
quark mass to be 171.8±7.2 GeV/c 2 . 



Appendix C 

DETAILS OF THE 
ALL-HADRONIC AND 
DILEPTON MASS ANALYSES 

C.l All-hadronic channel 

A reevaluation of the systematic uncertainties on the measurement 
of the top quark mass in the all-hadronic channel as reported in Ref. |ID[ 
has shown that some of those estimates were overly conservative. Since that 
publication further studies of the systematic uncertainties have led to better 
procedures, which we now apply to all channels. The systematic uncertain- 
ties which have been revised include: initial and final state radiation, fitting 
procedure, and jet energy scale. These revisions are discussed below. 

The contribution due to uncertainty in modeling initial and final 
state hard radiation was 8.0 GeV/c 2 . To evaluate this uncertainty, standard 
HERWIG ti events were compared to samples which were constructed to have 
smaller and larger fractions of events in which one or more of the final state 
jets did not match any of the daughter quarks from the ti decay. The most 
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evident difference between the samples was that the width of the reconstructed 
mass distribution broadened as this fraction increased. On the other hand, 
simulated experiments showed only a very small shift in the fitted top quark 
mass. The systematic uncertainty was evaluated as follows. We generated 
two samples of simulated experiments, (a) one using the default HERWIG 
templates, and (b) one using templates which were constructed to have 90% 
of events containing one or more jets that were not matched to the daughter 
quarks from the ti decay. For the default HERWIG sample, approximately 
60% of events have one or more jets not matched to a quark from the ti 
decay. In both cases, we evaluated the median and the rms width of fitted top 
quark masses from the simulated experiments. The systematic uncertainty 
was taken to be the quadrature difference of the widths between samples (a) 
and (b). This number was then added in quadrature with the small shift in 
the median mass which was observed between sample (a) and (b). Essentially 
all of the 8.0 GeV/c 2 uncertainty was from the increase in the width of the 
distribution of sample (b). Further studies show that the change in width of 
the reconstructed mass distribution with increased radiation is reflected in the 
statistical uncertainties returned by the fits for simulated experiments; thus 
the statistical uncertainty obtained from our fitting procedure for the data 
sample already takes into account this effect. A reevaluation, using the same 
procedure as described in Section p.2|, results in a contribution from this source 



of 1.8 GeV/c 2 |62 



Another large source of systematic uncertainty (5.2 GeV/c 2 ) came 
from the effect of selecting the second-best rather than the best kinematic fit 
to each event. A smaller contribution came from considering different ways 
of interpolating between likelihood values at discrete top mass values in order 
to find the maximum likelihood point. A third contribution came from the 
finite Monte Carlo statistics that provided the expected reconstructed mass 
distributions at different top mass values. The first two contributions are no 
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longer identified as sources of significant systematic uncertainty since they 
concern the robustness of the chosen method. The contribution from Monte 
Carlo statistics, of 0.3 GeV/c 2 , remains. 

The jet energy scale uncertainty was determined to be 5.4 GeV/c 2 . 
Part of that (3.7 GeV/c 2 ), was due to differences in the calorimeter energy 
scale between two versions of the detector simulation. The source of this 
uncertainty was later corrected. As a result, the 3.7 GeV/c 2 contribution to 
the uncertainty was eliminated. 

A small reorganization of the contributions has occurred, which we 
mention in order to avoid any confusion in a comparison with Ref. [L0||. The 
soft gluon uncertainty (3.0 GeV/c 2 ) has been moved from the "gluon radiation 
and fragmentation effects" to the "jet energy scale" category. The Monte Carlo 
generator uncertainty (0.8 GeV/c 2 ) has been assigned its own category. The 
result is a new systematic uncertainty of 5.7 GeV/c 2 , with a breakdown into 
different contributions as listed in Table |10.1| . 



C.2 Dilepton channel 

The top quark mass measurement in the dilepton channel uses eight 
observed events that pass the standard selection criteria used for the dilepton 



channel |12], B2|. The criteria require that the leptons have opposite charges, 
that there be at least two jets per event, and include cuts on the missing 
transverse energy and the lepton transverse energies. 

This measurement involves two steps: a top mass estimate is ob- 
tained for each event, and then a likelihood fit, which allows for the presence 
of background, gives an overall best estimate of the top quark mass. The 
second step is similar to that in the lepton+jets topology, but the first step is 
appreciably different. 

In order to get a mass estimate for an individual event, we determine 
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a weight distribution as a function of an assumed top mass, m t . First, we 
assume that the event originates from tt production and decay, that the leading 
two jets are b jets from top decay, and that the leptons (e or ji) are from 
associated W^-boson decays. Next, we assume a value for the top mass, m t , 
assume pseudorapidity values, i]i and r/2, for the two neutrinos, and solve for 
the two neutrino momenta. In general there are eight solutions because of a 
quadratic ambiguity in each neutrino's longitudinal momentum and a choice 
of pairing leptons with jets. For each solution, we denote as FJ T P the vector 
sum of the solution's neutrino transverse momenta. Then we assign a weight 
to each solution according to how well FJ T P agrees with the event's measured 
missing transverse energy, F^ T "\ as follows: 

g(m t , Vi, m) = exp (- 2a 2 j x ex P I \a 2 J ( ^ 

where a is the resolution in each component (x and y) of the measured un- 
clustered transverse energy (see below). The experimental resolution in jets 
and leptons is taken into account by sampling the measured quantities many 
times according to their resolutions. That is, for each set of assumed m t , r/i 
and 7/ 2 values a weight is calculated many times, and the sum is accumulated. 
For each assumed mt value, 100 pairs of r/i and r/2 values are assumed in turn, 
and the summed weights are again summed, to give a final summed weight, 
f(mt), at any m t value. The r/i and r/2 values are drawn independently from 
a Gaussian distribution with unit width and centered at 0.0, as predicted by 
HERWIG Monte Carlo simulations. Thus all the uncertainties on the F^ T mea- 
surement are taken into account, except for the resolution of the unclustered 
energy. We use a = A-Jn GeV, where n is the number of interactions in the 
event and comes from studies of low-luminosity minimum-bias events. 

For each event, m t values in the range 90 to 290 GeV/c 2 , in 2.5 
GeV/c 2 steps, were assumed in order to give a f(m t ) distribution. This distri- 
bution is used to determine a top mass estimate, as follows. The position of 
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the maximum value, f(m t ) max , is denoted by M max . The first points on either 
side of M max that have f{m t ) < f(jnt) m ax/2 are denoted by Mi and M 2 . The 
average of Mi and M 2 is taken as the top mass estimate. 

The f{m t ) distributions, normalized to unity, for the eight events are 
shown in Fig. |C.1| . The eight events, with their lepton identifications, numbers 
of jets (with uncorrected transverse energy greater than 10 GeV and pseudo- 
rapidity in the range —2.0 to +2.0 units), and estimated top masses are given 
in Table 



Run Event leptons Nj et Top mass \og(P ev ) 



41540 


127085 


e" 


> + 


2 


158.8 


0.47 


45047 


104393 


e H 


-fi 


2 


180.0 


1.82 


47122 


38382 


e H 


>- 


2 


176.3 


1.40 


57621 


45230 


e H 


>" 


2 


156.3 


2.20 


66046 


380045 




>~ 


4 


172.5 


-5.20 


67581 


129896 


e H 


>~ 


2 


143.8 


0.44 


68185 


174611 


e" 


^e- 


2 


161.3 


4.10 


69808 


639398 


e" 




3 


170.0 


3.50 



Table C.l: Information on the eight candidate dilepton events used in the 
dilepton mass analysis. Shown are: the run and event numbers, the types 
of leptons in each event, the number of reconstructed jets (with uncorrected 
Pt > 10 GeV/c and \r)\ < 2), and the top mass estimates for each event. Also 
listed is \og(P ev ), where P ev is the sum of all the weights for the event divided 
by the number of resolution samplings used. 

It is useful to define a variable, P ev , as the sum of f(m t ) over all 
assumed m t values, divided by the number of resolution samplings used. The 
latter number is 1500 for data and 30 for Monte Carlo events. This variable 
gives an indication of how easily an event can be fit to the tt decay hypothesis. 
The log(P e „) distribution of simulated it plus background events is shown in 
Figure |C.2| . The tt events are from the HERWIG simulation with a top quark 
mass of 175 GeV/c 2 . The log(P e „) values for the eight observed events are 
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listed in Table |C1| and are indicated by arrows in Figure |C2| . The data points 
all lie within the range spanned by the simulated distribution. In the simulated 
events, 0.7% have log(P e „) < —5.2, the value for the lowest data point, so the 
probability for an eight-event sample to have at least one event at —5.2 or 
lower is 5%. 

In Ref. it was noted that the same method could be applied to 
events in the lepton+jets topology that had two SVX-tagged jets. In such 
events the two untagged jets (of the four highest Et jets) are assumed to 
result from IV-boson decay, and in order to mimic a IV-boson leptonic decay 
one of those jets is treated as a lepton (electron or muon) and the other as 
a neutrino. In the following we took the jet with lower Et as an unobserved 
neutrino and recalculated FJ T m for the event. Then the above dilepton method 
was applied. 

The five events in the SVX Double sample were fit with this method. 
A top quark mass value of 181.5 ± 12.6 GeV/c 2 was obtained. This value 
has to be compared with the value shown in Table [D] of 170.0tg;g GeV/c 2 , 
a difference of 11.5 GeV/c 2 . In order to understand the difference between 
the two methods a comparison was made in a Monte Carlo study that used 
a sample of approximately 1300 simulated lepton+jets ti events with M top = 
175 GeV/c 2 and with two jets having SVX tags. The distribution of the 
reconstructed mass from the standard lepton+jets kinematic fit is shown in 
Figure C.3| (a). Also shown is the top mass estimate per event with the pseudo- 
dilepton method described above. The two distributions are similar. The 
medians are 170.5 GeV/c 2 and 170.9 GeV/c 2 , and the widths are 21.4 GeV/c 2 
and 23.4 GeV/c 2 , respectively for the kinematic fit and the dilepton methods. 
Here the widths are one-half the separation of the 16 th and 84 th percentiles 
in the distributions. As expected, the dilepton method gives a slightly wider 
distribution. In Figure |G3j (b) the mass difference between the two methods 
is plotted for each event. The width of this distribution is 24.3 GeV/c 2 . This 
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shows that the shift of 11.5 GeV/c 2 found for the five SVX Double events 
using the two methods is well within expectation. 

This study shows that fitting the dilepton events, which are undercon- 
strained, using the technique described here is just as valid and precise as the 
completely constrained 2-C fit used for the lepton+jets sample. In addition, 
if we calculate the statistical correlation between the two methods, we obtain 
a correlation coefficient of 0.36, i.e., fitting the SVX Double events with this 
technique could improve the statistical uncertainty on the mass determination 
from this channel. 
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Figure C.I: Weight distribution f(m t ), normalized to unity, for the eight 
observed dilepton events. 
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Figure C.2: Predicted distribution of log(P e „), the total weight sum per reso- 
lution sampling, for the expected ti and background event mix in the dilepton 
sample. The arrows indicate the values for the eight observed events. 
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Figure C.3: (a) Top mass estimates for simulated tt events (top mass 175 
GeV/c 2 ) with two SVX tags using the lepton+jets kinematic fit (dashed) and 
the pseudo-dilepton (solid) methods, (b) The difference per event between the 
top mass estimates from the two methods. 
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